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ABSTRACT 

We present a detailed analysis of the XMM-Newton and Chandra observations of Abell 2626 focused 
on the X-ray and radio interactions. Within the region of the radio mini-halo (~ 70 kpc), there are 
substructures which are probably produced by the central radio source and the cooling core. We find 
that there is no obvious correlation between the radio bars and the X-ray image. The morphology 
of Abell 2626 is more complex than that of the standard X-ray radio bubbles seen in other cool 
core clusters. Thus, Abell 2626 provides a challenge to models for the cooling flow - radio source 
interaction. We identified two soft X-ray (0.3-2 keV) peaks with the two central cD nuclei; one of 
them has an associated hard X-ray (2-10 keV) point source. We suggest that the two symmetric radio 
bars can be explained by two precessing jets ejected from an AGN. Beyond the central regions, we find 
two extended X-ray sources to the southwest and northeast of the cluster center which are apparently 
associated with merging subclusters. The main Abell 2626 cluster and these two subclusters are 
extended along the direction of the Perseus-Pegasus supercluster, and we suggest that Abell 2626 is 
preferentially accreting subclusters and groups from this large-scale structure filament. We also find 
an extended X-ray source associated with the cluster SO galaxy IC 5337; the morphology of this source 
suggests that it is infalling from the west, and is not associated with the southwest subcluster, as had 
been previously suggested. 

Subject headings: cooling flows — galaxies: clusters: general — galaxies: clusters: individual 
(Abell 2626) — intergalactic medium — radio continuum: galaxies — X-rays: 
galaxies: clusters 
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1. INTRODUCTION 

The central regions of clusters of galaxies are among 
the most interesting and physically active areas in the 
Universe. In many clusters of galaxies, there is a central 
peak in X-ray surface brightness, with a central cool- 
ing time which i s less than the age of the cluster (e.g., 
IChen et al.ll2007[) . In the classical cooling flow picture, 
the rapid cooling of gas in this region leads to the loss of 
central pressure (|Fabianl H994). If there is no additional 
heating, the cooled gas would flow into the center in a 
cooling flow. The gas should continue to cool to very 
low temperatures. However, recent high resolution X- 
ray spectra from XMM-Newton and Chandra show that 
there is not enough cooler X-ray gas to be con s istent 
with a classical cooling fl ow (|Peterson et al.ll200ll l2003t 
iPeterson fc Fabianl[200a ). 

The lack of sufficient cooler materials probably means 
that there is an additional heating source within the cen- 
tral cooling core region. Very often, the clusters host 
central dominant galaxies which are also strong radio 
sources. In some of the clusters, there is strong evidence 
that the central radio sources are interacting with the 
surrounding X-ray emitting plasma. The radio sources 
blow "bubbles" in the X-ray emitting gas, displacing the 
gas and compressing it into shells around the radio lobes. 
At the same time, the radio sources are confined by 
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the X -ray emitting gas (tFabian et al.f 2000: Blanto n et alJ 
2001). In some clusters, "ghost bubbles" at larger radii 
are seen whi ch are weak in radio emission except at low 
frequencies (jMcNamara et al.l l2001t IClarke et al J 12005( 1 . 
They are believed to be previous eruptions of the radio 
sources, which are episodic. In some cases, buoyantly ris- 
ing bubbles may entr ain cooler X-ray gas from the centers 
of the cool ing cores (jChurazov et alJl2001t iFuiita et al.1 
l2002t 120041 ). Some radio sources previously classified as 
cluster merger radio relics may actually be displaced ra- 
dio bubbles from the central radio sources (jFuiita et al.1 
2002] 120041 ). In some recent studies, there is even evi- 
dence of radio jets restarted at a different position angle 
from a previous outburst, w hich may be due to a precess- 
ing supermassive blackhole dGitti et al.ll2006h. Wi th high 
resolution 3D simulations, IVernaleo fc Reynolds! (|2006l ) 
have shown that multiple bursts from the AGN along 
the same direction fail to produce a long-term balance 
of heating and cooling. The reason is that a low-density 
channel is evacuated by a previous jet, while the subse- 
quent jets can pass through the channel freely without 
heating the cooling core region. It has been suggested 
that jet precession is a possible solution to this problem. 
Alternatively, turbulence, rotation, or bulk motions in 
the intracluster medium (ICM) might also move or dis- 
rupt the channel made by the radio source, and keep 
the energy from j ust going down the same channel (e.g., 

iHeinz et all 120061 ). 

Abell 2626 (z = 0.0573. lAbelTeniI]fl98l hosts a mod- 
erate cooling flow of ~ 50 M yr _1 based on Einstein 
data (|White et alJll997f ). A double nuclei cD galaxy is 
sitting at the center of the clust er. From high resolution 
VLA images, iGitti et al.l (J200l recognized that there is 



2 



Wong et al. 



an unresolved radio core coincident with the central cD 
galaxy. There is a small jet-like feature extended towards 
the southwest direction from the unresolved radio core. 
In addition, there are two unusual amorphous, symmet- 
ric radio "bars" running parallel at opposite sides of the 
central cD galaxy. These compact features are distinct 
from and emb edded in dif f use ex tended radio emission 
(a mini-halo). iGitti et "all (|2004f ) believed that the ra- 
dio bars may represent an earlier evolutionary stage of 
jets injected by the central source. Thus, Abell 2626 
is a good candidate for studying the interaction be- 
tween the X-ray emit ting intracluster medium and the ra- 
dio emitting plasma (iR izza et al.l l2000t IGitti et al.lf2004 
iMarkovic etaill2004D . Previous ROSAT X-ray and VL A 
radio observations suggested that there is an X-ray excess 
which is spatially correlated with the radio source, but 
failed to find stron g X-ray defici t s (hol es) in Abell 2626 
(|Rizza et al.ll2000ft . iRizza et all (|2000l ) carried out nu- 
merical simulations, and showed that even a weak jet 
(internal Mach number = 5) should produce a significant 
hole along the line-of-sight to the radio jet. They sug- 
gested that the lack of X-ray holes may be due to mixing 
of thermal and radio plasma in the region, although their 
numerical simulations did not distinguish between these 
two components. They also ran simulations in which the 
jet turned off, which resembled the observations of weak 
X-ray holes. Another possibility they noted is the impre- 
cise model they used to subtract out from the symmetric 
X-ray emission to obtain a residual image. They searched 
for X-ray and radio correlations from the residual image 
but failed to find strong X-ray deficits (holes) displaced 
by the radio plasma. With the excellent spatial resolu- 
tion provided by Chandra and deeper observations pro- 
vided by both Chandra and XMM-Newton, we are able 
to study the complicated central region of Abell 2626 in 
more detail. 

While many radio sources are associated with AGNs in 
cluster galaxies, there exist other types of extended radio 
emission associated with the diffuse intracluster medium 
without clear galaxy hosts. These extended radio sources 
can b e divided into cluster-wide halos, relics, and mini- 
halos (|Feretti k. Giovan nini 1996 .120071) Abel l 2626 con- 
tains a candidate mini-halo (jGitti et al.l l2004). Cluster- 
wide halos and relics are large in scale (~ 1 Mpc), with 
halos located in the centers of clusters, while the relics 
are generally located at the periphery. Very often, these 
two types of large scale diffuse radio emission are found in 
clusters without cooling cores. This indicates a possible 
origin in cluster mergers, which are believed to disrupt 
coolin g cores and to re-accelerate rad io-emitting parti- 
cles (|B rune M 12004 iBlasi et al.1 120071 ) . However, some 
cooling core clusters with dominant radio galaxies host 
mini-halos at the cluster center with extended radio emis- 
sion up to 0.5 Mpc. The morphology of the mini-halos 
is different from that of the radio bubbles at the cen- 
ters of many cool core clusters, which show a strong X- 
ray/radio interaction in which the radio lobes displace 
X-ray gas. This hints at a different transport mecha- 
nism (not jets) to explain the extended property of the 
mini-halos. The typical time scales for radiative losses 
by relativistic electrons (~ 10 7 — 10 8 yr) in cooling cores 
with strong magnetic fields are much shorter than the 
diffusion times (~ 10 9 yr) for the electrons to be trans- 
ported from the center of the cooling core out to a frac- 



tion of a Mpc. Hence, the diffuse radio emission from 
mini-halos suggests that an in-situ re-acceleration mech- 
anism is needed for the electrons. The origin of the 
mini-halos has remained unclear, but it has been sug- 
gested that the interaction of the central radio galaxy 
and the intracluster medium (r e-) accelerates part icles to 
relativistic energies. Recently, IGitti et all (|2004D mod- 
eled Abell 2626 with an electron re-acceleration model 
driven by MHD turbulence amplified by the compression 
of the magnetic field in the cooling core region. They 
have shown that the model agrees with various obser- 
vational constraints such as the observed radio bright- 
ness profile, the integrated synchrotron spectrum, and 
the radial steepening of the radio spectrum. They con- 
clude that only a tiny fraction (~ 0.7%) of the maxi- 
mum power that could be extracted from the cooling core 
(based on the standard cooling flow model) is needed for 
re-acceleration. Since the re-acceleration model for the 
origin of the mini-halo assumes a physical connection be- 
tween thermal and relativistic plasma, it is essential to 
perform new X-ray observation to derive accurately the 
properties of the ICM (e.g., cooling radius, mass accre- 
tion rate, as well as morphology of the X-ray emission). 
On another hand, iPfrommer fc Enfilinl (|2004l ) suggested 
that the re-acceleration is of hadronic origin which should 
be best tested by 7-ray observations. 

In this paper, we present a detailed study of Abell 2626 
from XMM-Newton and Chandra observations focused on 
the X-ray and radio interaction. The paper is organized 
as follows. In § [21 we describe observations and data re- 
duction. The X-ray images are analyzed in § [3J with an 
emphasis on the X-ray morphology of the overall region, 
the central region, and the SO galaxy IC 5337 in § 13.11 In 
particular, two extended X-ray emission regions possibly 
associated with subclusters are discussed in § 13.1.11 The 
point sources detected in the X-ray data are discussed in 
§ 13.21 The X-ray surface brightness profile is studied in 
§ 13.31 Various substructures identified on residual maps 
are presented in § 13.41 Spectral analysis is presented 
in § [H In particular, cluster profiles, mass deposition 
rate, cooling and thermal conduction time scales, hard- 
ness ratio maps, possible AGNs associated with the cD 
galaxy IC 5338 and the SO galaxy IC 5337 are discussed in 
this section. Throughout the data analysis sections (§§[3] 
& S]), some discussions related to the interpretation of 
the results are presented to help readers understand the 
data. We then discuss the kinematics of the SO galaxy 
and the cooling flow interaction with the central radio 
source in §§[5]&[6]in detail, respectively. We summarize 
our work in § [7] We assume Hq = 71 km s _1 Mpc -1 , 
£1m = 0.27 and £l\ — 0.73 throughout the paper. At 
a redshift z = 0.0573, the luminosity distance of Abell 
2626 is 252.8 Mpc, and 1" = 1 10 kpc. The virial ra- 
dius is approximately 1.53 Mpc (M arkevitch et al.lll998h 
for the average cl uster temperature of Tx — 3.0 keV 
(|White et al.lll99l . 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. XMM-Newton Data 

Abell 2626 was observed (Obs. ID 0148310101) on 2002 
December 28 by XMM-Newton for ~ 38 and ~ 42 ks with 
the EPIC PN and MOS cameras, respectively. Extended 
Full Frame mode (PN) and Full Frame mode (MOS) were 
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used with thin optical filters. XMM-Newton data reduc- 
tion was done with SAS version 5.4.1. All the data were 
reprocessed using the tasks epchain and emchain. We 
included only events with FLAG=0 in our analysis. We 
chose only events with PATTERN = 0-4 (0-12) for PN 
(MOS) data. 

Background flares were rejected by performing a 2a 
clipping of the 100 s binned light curves in the 12-14 
keV and 10-12 keV band for the PN and MOS cameras, 
respectively, which yielded the out-of-time corrected ex- 
posure times of 35,563 s, 40,300 s and 40,600 s of cleaned 
data for the PN, MOS1 and MOS2 cameras, respectively. 
Blank-sky background files collected by D. Lumb 4 were 
used during the analysis. The observed background lev- 
els, estimated from the total count rate in the whole 
field of view in the hard band (10 - 12 keV and 12 - 
14 keV for MOS and PN, respectively), are 1.07, 0.99 
and 1.26 times those in the blank-sky backgrounds for 
MOS1, MOS2 and PN, respectively. We have renormal- 
ized the blank-sky backgrounds accordingly. The uncer- 
tainty of the PN background is larger than that for the 
MOS detectors. The effect of the larger PN background 
uncertainty and some other background uncertainties on 
the spectral analysis are addressed in § 14.1.11 The im- 
age and hardness ratio map analyses were done using the 
MOS data only (Stil3lfcl43i and we checked that the re- 
sults of the analysis would not be affected by a change in 
the background normalization of ±10%. The vignetting 
effect for the image and hardness ratio map analysis has 
been corrected using the exposure maps. 

2.2. Chandra Data 

A Chandra observation (Obs. ID 3192) of ~ 25 ks was 
made on 2003 January 22 in Very Faint mode. The center 
of the cluster was positioned 3' from the edge of the S3 
chip (midway between the node boundaries). The focal 
plane temperature was —120 C. 

The Chandra data reduction was done with CIAO ver- 
sion 3.1, using calibration products from CALDB 2.28. 
The latest gain file (CALDB GAIN 2.28) and geometry 
file (CALDB 2.9) at the time of reduction were used. The 
reduction was done in Very Faint mode to remove more 
particle background. Background flares were rejected by 
performing the standard clipping using the lc_clean rou- 
tine written by M. Markevitch 5 . This yielded 22,959 s of 
cleaned data. Since the whole field-of-view is dominated 
by the cluster X-ray emission, internal background can- 
not be used. The blank-sky background files collected 
by M. Markevitch were used during the analysis. The 
difference between the observed and the blank-sky high 
energy (9-12 keV) background count rates was less than 
1%, and this small difference was removed by renormaliz- 
ing the background. The influence of the background on 
the spectral fitted measurements (temperature profile) 
has been checked by varying the background normal- 
izations by ±10%, and the results are not significantly 
changed. To check if there might be a significant dif- 
ference between the Galactic background in the source 
and blank-sky observations, another thermal component 
with a fixed low temperature of 0.2 keV ha s been added 
(jMarkevitch et al.ll2003l : iLumb et aT1l2002f ). Normaliza- 

4 ftp: // xmm. vilspa. esa.es/pub/ ccf / constituents /extras /background / 

5 http:/ /cxc. harvard.edu/contrib/maxim/acisbg/ 



tion of the additional component has been allowed to be 
negative for the case of more Galactic emission in the 
blank-sky observations compared to the cluster observa- 
tion. All the results are basically unchanged. Fixing 
the temperature of the additional component to be the 
outermost temperature that can be determined from the 
Chandra field of view also gives essentially the same re- 
sults. We conclude that the Chandra results are insensi- 
tive to uncertainties in the background. 

3. X-RAY IMAGE ANALYSIS 
3.1. X-ray Morphology 
3.1.1. Global morphology 

The background subtracted, exposure-corrected, adap- 
tively smoothed mosaic of the XMM-Newton EPIC 
MOS1 and MOS2 images is shown in Figure □ The 
PN image was not included because of its different spec- 
tral response and because the chip gaps and bad columns 
produced cosmetic artifacts in the smoothed image. The 
image was smoothed to a signal-to-noise ratio of 3 per 
smoothing beam. 

The global X-ray image of Abell 2626 appears to be 
roughly azimuthally symmetric. The brightest region at 
the center corresponds to the position of the cD galaxy 
IC 5338. There is a source associated with the cluster 
SO galaxy IC 5337 which is 1.'3 west of the center of 
the cluster. Excess X-ray emission is seen to the south- 
west of the cluster, exten ding from ~1' out to ~7' (see 
§ELH). iMohr et all (|1996h identified a subcluster to the 
southwest of this region based on optical observations 
of the galaxy population. The subcluster identified is 
centered roughly at the southwest edge of the XMM- 
Newton field of view, but the distribution of the galaxies 
covers a reg i on wh ere the excess X-ray emission is seen. 
IMohr et~al . ( 1996) show that the subcluster is bound to 
and falling into the main cluster. Our XMM-Newton im- 
age shows, for the first time, extended X-ray emission 
possibly associated with the subcluster in the southwest 
direction. The extended X-ray emission can be viewed on 
a background subtracted, exposure corrected, adaptively 
smoothed image. The location of the extended emission 
possibly associated with the subcluster is circled at the 
southwest corner in Figure [TJ The extended emission at 
the southwest corner is not obvious, but the slight devia- 
tion from spherical symmetry in that region can be noted. 
The shape of the extended emission region is difficult to 
characterize due to the low signal-to-noise ratio near the 
edge of the cluster. By inspecting various images, we 
chose a circular region centered on the highest excess X- 
ray emission, with a radius chosen to be small enough so 
that the region did not overlap with the nearly circular 
cluster emission of Abell 2626. Our choice for this region 
is justified by the hardness ratio map below (Figure [14)) 
which shows that the extended region chosen appears to 
be harder than its surrounding. Compared to the X-ray 
emission at the same radius from the main cluster cen- 
ter (the cD galaxy IC 5338) with all point sources and 
other extended emission excluded, the excess emission is 
significant at the 4.6 a level. Our residual images (§ 13.41 
below) also indicate the possible signatures of a merger. 

Another extended X-ray source can be seen about 7' 
northeast of the center of Abell 2626 (Figure [T]) . The 
criteria used to define the northeast region are the same 
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as the southwest one. This feature is significant at the 
11.6 a level. It should be noted that there is a chip gap 
at about 5.'5 from the center which makes the extended 
X-ray source appear to be more distinct than it actu- 
ally is. No cluster or group structure has been identi- 
fied previously in this region. However, Abell 2626 is 
known to be associated with the Perseus-Pegasus su- 
percluster, a filament of clu sters of galaxies e xtending 
for as much as ~ 300 Mpc dBatuski , fc Burns! ll985allbT 
IZucca et all Il993t lEinasto et all 120011) . The southwest 
to northeast extension to the structure in our X-ray im- 
age around Abell 2626 is elongated along the Perseus- 
Pegasus filament. This may indicate that the cluster is 
preferentially accreting subclusters and groups from this 
supercluster. 

3.1.2. Central Region 

The background subtracted, exposure corrected, adap- 
tively smoothed Chandra image of the central 3.'5x2f4 is 
shown in Figure [2] (left panel) . The image was smoothed 
to a signal-to-noise ratio of 3. The brightest region of the 
X-ray emission in the cluster is centered on the southwest 
nucleus of the central cD galaxy. This source is possibly 
extended in both soft (0.3-2 keV) and full band (0.3-10 
keV) images when detected by wavedetect in CIAO, but 
nothing was found by wavedetect in the hard band (2-10 
keV). However, examination of the raw soft band (0.3-2 
keV) image shows two peaks superposed on the two op- 
tical nuclei of the central cD galaxy (lower left panel of 
Figure [3]). There is also a hard point source seen in the 
raw hard band (2-10 keV) image (lower right panel of 
Figure [3]) . The position of this hard point source agrees 
with the southwest cD nucleus. The hard point source is 
separated from the peak of the radio core by less than 1". 
This source can also be seen in the hardness ratio map 
(Figure [TBI below) . The hard point source and its likely 
association with the AGN at the center of the southwest 
cD nucleus and radio core are discussed further in § 14.71 
This source is located along an arc of X-ray emission ex- 
tending from the northeast to the west to the south of 
the cluster center. 

Previous ROSAT X-ray and VLA radio observations 
suggest that there is enhanced X -ray emission spat ially 
coincident with the radio source (R 

izza etaDEOOS). In 
Figure [2] (left panel), the contours from th e VLA 1.5 
GHz B-array radio image are shown in green (|Gitti et al.l 
2004). The position of the radio core agrees with the 
brightest region of X-ray emission. In particular, the 
radio core is centered at the southwest cD nucleus rather 
than the northeast one. There is an arc of X-ray emission 
which may correspond to the curved radio jet feature to 
the south of the core. 

However, the most unusual features of the radio image 
are the two elongated radio bars to the north and south 
of the radio core. We find there is no obvious correla- 
tion between the two radio bars and the Chandra X-ray 
image in Figure [2] (left panel). No obvious X-ray deficit 
(holes) is found on the Chand ra X-ray image, wh ich is 
consistent with the finding of iRizza et al.1 |2000). We 
have generated residual maps to further investigate our 
finding, and the results will be presented in § 13.41 Our 
interpretation is that the radio bars are thin tubes (see 
§0 below). 

One interesting feature is that there is an X-ray ex- 



cess "tongue" which extends from the central cD galaxy 
to the southern radio bar, and there is a similar but 
weaker tongue in the northern direction as well (upper 
right panel of Figure [3]) . The tongues are best seen in 
the residual map in Figure [8] of § 13.41 

Figure [2] (left panel) also show s radio contours ( white) 
from the VLA 1.5 GHz C-array (jGitti et all 120041 ) over- 
laid on the Chandra image showing the radio mini-halo 
in Abell 2626. The mini-halo is roughly confined within 
60" with a diamond shape. The distance from the cen- 
ter to the NE or SW corners of the diamond is about 
50", while to the NW or SE corners it is about 70". A 
discussion on the mini-halo will be presented in § [6] 

3.1.3. The SO Radio Galaxy IC 5337 

About west of the central cD galaxy, there is an 
X-ray source associated with the SO galaxy IC 5337 (left 
panel of Figure [2]) . This source is clearly extended in the 
soft band image, but a hard band image shows that there 
is a point source coincident with the nucleus of IC 5337. 
The possibility that the SO nucleus is an AGN will be 
discussed below f§ l4.7[) . To the south of IC 5337, there is 
another X-ray point source which is not associated with 
any known galaxy. As indicated in Figure [2J there is a 
nearby radio source which may be associated with this 
X-ray source. 

The right panel of Figure [2] shows the bow-shock-like 
shape of the extended X-ray emission associated with the 
SO galaxy IC 5337. The image is created by removing and 
replacing the southern point source by the average sur- 
rounding X-ray intensity, and then adaptively smoothing 
to a signal-to-noise ratio of 4. 

The SO galaxy IC 5337 also is a radio source (Figure[2]). 
Radio contours (jGitti et al.ll2004f ) are overlaid in Figure[2] 
(right panel). There is a component centered on the SO 
galaxy, and a component located at a position of about 
190° (measured from north to east) of the SO galaxy on 
the 1.5 GHz B-array map. The 1.5 GHz C-array map 
shows three radio tails in the south, southwest and west 
directions. The kinematics of the SO galaxy will be dis- 
cussed in § [5] 

3.2. Source Detection 

We used the XMM SAS task eboxdetect to detect a 
total of 90 distinct point sources on the XMM images; 
the detections were done separately on each of the three 
detectors, and then the source lists were combined. The 
sources were confirmed by inspection, and some low-level 
detections at the edge of the FOV or bad pixel gaps were 
removed. Here we discuss only the sources with possible 
optical counterparts within 20". In fact all optical coun- 
terparts identified are within 11". We adopt the NED 
source with an accurate position which is closest to the 
X-ray position, The possible associations of X-ray point 
sources with AGNs in the central cD galaxy IC 5338 
and the SO galaxy IC 5337 were noted above 13.1.21 &; 
I3.1.3p . and will be discussed in more detail below (§ 14.7 ). 
Nine X-ray sources are included in the list (Table [T|. 
Figure [4] is an optical image from the Digital Sky Sur- 
vey with the 9 X-ray sources marked. Two of these are 
the central cD galaxy IC 5338 and the SO galaxy IC 5337 
discussed above. Interestingly, three of the X-ray sources 
correspond to galaxies which are located in the region of 
extended X-ray emission about 7' northeast of the center 
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of Abell 2626 (Figure [T]). However, one of these possible 
optical IDs (2MASX J23365722+2114032) has a redshift 
of 0.038, and thus is not directly associated with Abell 
2626. 

On the Chandra image, 14 point sources were detected 
by wavedetect in CIAO using the default settings. Of 
these, only the central cD IC 5338 galaxy and the SO 
galaxy IC 5337 have possible identifications in NED with 
positions that agree with the X-ray to within 10". Ta- 
ble lists the positions of the 14 Chandra X-ray point 
sources. Figure [5] is an optical image from the Digital 
Sky Survey with the locations of the 14 X-ray sources 
marked. 

Unless specified, all the sources detected are excluded 
in the analysis of extended emissions in the subsequent 
sections. 

3.3. Surface Brightness Profile 

The azimuthally- averaged X-ray surface brightness 
profile was extracted separately for the Chandra and 
XMM-Newton data. Except for the two extended X-ray 
clumps to the NE and SW and the structure in the very 
center, Abell 2626 appears to be a nearly relaxed cluster. 
The profiles were fitted using a double beta model profile 
of the form 



S(r) = Si(0) 



1+ ( — 

' ri 



-3ft+l/2 



(1) 



which has been shown to provide a good fit to relaxed 
clusters with cooling cores (|Xue fc Wdl2000h . 

The background subtracted and exposure corrected 
surface brightness profile from the XMM-Newton MOS1 
and MOS2 data is shown in Figure[6] The PN data were 
not included since the spectral response is different which 
complicates the deprojection of the electron density pro- 
file below (§ 14. 3p . Also, the fitted surface brightness 
model was used to create residual images (§ 13.41 below), 
and including the PN image would produce cosmetic ar- 
tifacts due to the larger chip gaps and greater number of 
bad columns in the PN camera. All of the X-ray point 
sources except the central cD galaxy were excluded from 
the data. The extended X-ray emission around the SO 
galaxy was also excluded, as it produced a feature in 
the cluster surface brightness profile. We also excluded 
the extended X-ray emission regions to the NE and SW 
of the cluster which were described in and the criteria 
of defining the regions were explained in § 13.1.11 The 
annuli were chosen such that all the widths are larger 
than the FWHM of the PSF at that radius. With this 
choice, all the bins contain at least 1000 counts after 
background subtraction for the XMM-Newton MOS1 and 
MOS2 data. With these very good statistics, the error 
bars on the surface brightness are too small to be easily 
seen in the top panel of Figure [6l instead, we present 
the residuals in the lower panel. The parameters of the 
best-fit double beta model are Si(O) = 0.40t^? and 

S 2 (0) 
and r2 = 

02 = 0.731^;^. Visually, the double beta model fit seems 
adequate. However, with the very small errors on the 
measured surface brightness, the fit is actually not very 
good, with a reduced \ 2 of 2.3. In addition to the excess 



i.02 

n = 34.81^ 
128.5^2 arcsec, and ft = 0.91±g;|l and 



O.lOI^oi counts s 1 arcmin 2 , 



at the center, associated with the central cD galaxy and 
radio source, within 200" there are several other regions 
of excess or deficient emission compared to the model. 
These features are reproduced (albeit with poorer statis- 
tics) in the Chandra surface brightness data, and in the 
data from the XMM-Newton PN detector. These resid- 
uals may indicate that the central region is disturbed by 
effects of the central radio source beyond the current ex- 
tent of the source, or that the cluster is not completely 
relaxed due to a merger. 

3.4. Residual Images 

To better understand the origin of the residuals in 
the surface brightness fit and the dynamics of the in- 
tracluster gas, we created residual maps for both XMM- 
Newton and Chandra images (Figures [7] and |HJ). Sim- 
ulated azimuthally-symmetric X-ray images were pro- 
duced which exactly followed the double beta model 
fits to the surface brightness. These simulated images 
were subtracted from the adaptively smoothed images. 
Figure [7] shows the XMM-Newton and Chandra resid- 
ual images with the same field of view. Both residual 
maps show similar structures; the XMM-Newton image 
has great contrast due to the better statistics, while the 
Chandra image shows more detailed structures. Sub- 
structures can be identified with many of the residuals 
seen in the maps. A large region of excess emission can 
be seen about 110" southwest of the cD galaxy, which 
corresponds to the peak of the surface brightness resid- 
uals around that radius. Other regions of excess and 
deficit within 100" can also be identified from the resid- 
ual map of the Chandra image (right panel of Figure [7]). 
In the surface brightness residuals, the central excess cor- 
responds to the central bright point in Figure [7] (right 
panel). The trough around 15" corresponds to the dark 
ring around the central excess in the residual map which 
can be seen more clearly on FigureO The peak at around 
30" is probably due to the two bright regions due north 
and south in the residual map (labeled as two green cir- 
cles). The trough around 50" corresponds to the dark 
region indicated by a white polygon in the residual map. 
The peak around 110" corresponds to the far southwest 
excess labeled with a green polygon in the residual map. 

Most of the structures in the inner region occur on a 
scale similar to that of the radio sources in the center. 
In order to study the radio/X-ray interaction in more 
detail, Figure [8] shows radio contours from the VLA 1.5 
GHz C-array (white contours ) and the 1.5 GH z B-array 
radio images (green contours. IGitti et"al"1l2004h . overlaid 
on the central part of the Chandra residual map. In gen- 
eral, the diamond-shaped outline of the extended radio 
emission (the radio "mini- halo" ) corresponds to a region 
of reduced X-ray emission, and there are excess X-ray re- 
gions NNE and SSW of the edge of the mini-halo where 
the radio surface brightness is dropping rapidly. The 
mini-halo is roughly confined within 60". The distance 
from the center to the NE or SW corners of the dia- 
mond is about 50"; while to the NW or S E corners is 
about 70". Note that in the interpretation of lGitti et al.l 
(2004), the two radio bars are distinct from the diffuse 
radio emission, and in fact after their subtraction the 
morphology of the mini-halo be comes roughly circular 
(see Figure 6 oflGitti et al.ll2 004). However, in the resid- 
ual image in IGitti et al.l (|2004). we can still see some ex- 
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cess radio emission in the NW and SE corners. As noted 
previously, there is a significant excess associated with 
the radio core and the southern radio jet. The southern 
radio bar is a region of X-ray excess, but the northern 
bar shows an excess to the east and a deficit to the west 
in X-rays. The interpretation of the interaction between 
thermal and radio plasma will be discussed in § [5] in de- 
tail. 

The one major X-ray residual in the inner cluster which 
is not on the scale of the central radio emission is the 
extended excess to the SW. It is possible that this has 
been produced by the action of the central radio source. 
However, the SW excess is between the cD galaxy and 
the subcluster southwest of the main cluster center. It is 
possible that this excess is the result of this merger or a 
previous merger from this direction. There is no obvious 
sharp surface brightness discontinuity associated with 
this excess, which means that it either is not a shock or 
cold front, or that the merger is not occurring pr i marily 
in the plane of the sky. We note that iMohr et alj (1996) 
show that the line-of-sight velocity difference between the 
main cluster and the subcluster is ~ 2600 km s _1 , which 
is greater than the typical sound speed within clusters. 
If the excess is the result of a merger, this suggests that 
the merger is mainly along the line-of-sight. 

4. SPECTRAL ANALYSIS 

4.1. Temperature Profile 

The azimuthally averaged, projected temperature pro- 
file for Abell 2626 was extracted using both XMM- 
Newton and Chandra data (upper panel of Figure 
The spectra were fitted using an absorbed thermal model 
(WABS*MEKAL) with XSPEC Version 12.3.1. For the 
XMM-Newton data, all the spectra of the three cameras 
were fitted simultaneously. The energy bands used are 
0.6-10.0 and 0.4-7.0 keV for XMM-Newton and Chandra 
data, respectively. Using 0.6-7.0 keV and 0.5-10.0 keV 
bands for XMM-Newton and Chandra data gives essen- 
tially the same results. The spectra were grouped to have 
a minimum of 25 count s per bin. We fixed the ab sorption 
to be 4.2 x 10 20 cm -2 dDickev k, Lockmanlfl99fjh and the 
redshift to be 0.0573 (jAbell et al J 1 19891) . Freeing these 
values gives essentially the same results. Except the out- 
ermost data point where the reduced \ 2 value is ~ 1.6, 
all the other reduced % 2 values are < 1, which indicates 
that the fits are good. To ensure that the spectra are not 
affected by energy-dependent PSF effects, the annuli of 
the XMM-Newton data were chosen such that all of the 
widths are at least twice the PSF half energy width of 
the PN camera. 

The general trend is that the projected temperature 
increases from a central value of ~ 2.5 keV within 9" 
(~ 10 kpc) to a maximum of - 3.5 keV at 100" - 200" 
(110 — 220 kpc). Such a radial increase is normally seen in 
a cluster cooling core. The projected temperature then 
decreases to ~ 1.4 keV at ~ 540" (~ 600 kpc). Inter- 
estingly, there is a sharp drop in the projected temper- 
ature from ~ 2.9 keV at 210" (230 kpc) to - 2.0 keV 
at 270" (296 kpc) in the XMM-Newton data. Outside of 
this steep drop, the projected temperature gradient flat- 
tens. The radius of the drop corresponds to the distance 
from the cD galaxy to the outer edge of the large excess 
emission in the southwest direction. If the large excess 
emission indicates a merger signature, the temperature 



jump might indicate that gas is heated within that re- 
gion. Unfortunately, the jump cannot be confirmed with 
Chandra as the Chandra data does not extend far enough 
in all directions to reach the distance of the jump seen in 
the azimuthally averaged temperature profile. It only ex- 
tends far enough towards the northwest direction, where 
there are not enough counts to constrain the tempera- 
ture. 

In general, the XMM-Newton data give lower values 
for the projected temperature profile. There are at least 
three possible explanations for this difference. First, it 
might be due to the larger PSF of XMM-Newton. The 
softer photons from the brighter, cooler inner regions 
may be scattered to outer regions, hence making the 
outer regions appear cooler. This effect may be seen 
from the second bin outward, and should be more seri- 
ous in regions where the gradient of the surface bright - 
ness is greatest (i.e., near the center; [Markevitch|[2002). 
This might explain the discrepancy in the second and 
third bins of our XMM-Newton data. The second pos- 
sibility is that there might have been a steady and mild 
enhancement in the hard X-ray background during the 
Chandra observation which was not removed by the stan- 
dard cleaning procedure; this could make the Chandra 
emission appear hotter, particularly in outer, low surface 
bri ghtness regi o ns. S uch an effect has also been noted 
bv lMarkevitchl (|2002t ). We have tried to assess the im- 
portance of such a background enhancement by includ- 
ing a component for the hard backgro und, modeled as a 
power law with an exp onential cutoff (|Markevitchll2002l : 
iMarkevitch et all2003l ) in the fits of the Chandra spectra. 
However, this component did not significantly improve 
the fits, and the best-fit amplitude of the component 
was negligible. A third possibility is that the continuous 
accumulation of the contaminant layer on the Chandra 
ACIS optical blocking filter was underestimated. 

To study the true temperature of the cluster, we have 
also deprojected the temperature profile by using the 
XSPEC PROJCT model. For each annulus, the model 
sums up the appropriate fractions of emissions con- 
tributed from the outer annuli. The spectra of all the 
annuli were fitted simultaneously to get the deprojected 
temperature profile. The middle panel of Figure [9] shows 
the azimuthally averaged, deprojected temperature pro- 
file for Abell 2626 using both XMM-Newton and Chandra 
data. For the XMM-Newton data, all the spectra from 
the three cameras were fitted simultaneously. The XMM- 
Newton and Chandra data are consistent within the 90% 
confidence uncertainties, but the uncertainties of the 
Chandra data are much larger. Since the XMM-Newton 
data extends further out and has smaller error bars, we 
describe results obtained from the XMM-Newton in the 
following. We also adopt the XMM-Newton spectral fits 
(lower panel of Figure [9]) for the deprojection and the 
mass profile analysis below (§§ 14.31 l4~4l fc !4.5|) . 

The deprojected central temperature is slightly lower 
than the projected one because the latter is contaminated 
by the hotter emission along the line of sight. The differ- 
ence is small since the hottest region has a temperature 
only 1 keV higher than the central one. A more notice- 
able difference is that the temperature jump at around 
240" is more obvious in the deprojected profile. It drops 
from a temperature of ~ 3.5 keV to ~ 2.2 keV. The 
peak of the deprojected temperature profile is located 
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just within the temperature jump (~ 180 — 240"). 
4.1.1. Effects of XMM Background Uncertainties 

Since the cluster emission covers the whole field of view 
and the correct background cannot be estimated easily, 
the so called "double background subtrac tion" method, 
comm on in XMM-Newton data analysis (jArnaud et al.l 
2002), cannot be applied. The natural background sub- 
traction method is to perform a renormalized blank-sky 
background subtraction without the second step of ex- 
tracting a blank-field-subtracted spectrum in an emis- 
sion free region in the "double background subtraction" 
method. We use this "single background subtraction" 
method. Note that this method has also been applied to 
the study of a cluster with a similar redshif t to Abell 2626 
( Arnaud et al.l l2~001h . As mentioned in lArnaud et al.l 
(|2001h . our treatment does properly account for the 
cosmic-ray (CR) induced background. However, the con- 
tributions from variations in the local soft Galactic back- 
ground (which depends on the position on the sky) and 
the extragalactic X-ray background component (which 
depends on the absorbing hydrogen column density along 
the line of sight and cosmic variance) were not included. 
In our case, the emission from Abell 2626 usually dom- 
inates the local soft Galactic background, and the CR 
induced background usually dominates the hard X-ray 
extragalactic bac kground component. U sing the ROSAT 
All Sky Survey (|Snowden et alj Il997h . we have com- 
pared the soft X-ray emission (0.47-1.21 keV) around 
Abell 2626, and the exposure- weighted emission in the re- 
gions used to extract the blank-sky background, and the 
difference is less than 4%. The hydrogen column density 
towards Abell 2626 is about a facto r of two higher than 
that for most of the blank-sky fields (Dic key fc Lockmanl 
1990), but is not high enough to affect the cosmic X- 
ray background strongly. Also, as mentioned previously, 
the CR induced background usually dominates the hard 
X-ray extragalactic background component. Hence, the 
blank-sky backgrounds should represent the soft Galactic 
background and the hard X-ray background reasonably. 

We were concerned that, while the background renor- 
malizations for MOS1 and MOS2 are modest (< 7% 
correction), the renormalization of the PN background 
is larger (~ 26% correction). While the uncertainties 
in the background will affect most strongly the derived 
temperatures of the outer regions of the cluster, we em- 
phasize here that most of our focus in this paper is on 
the X-ray and radio interaction in the innermost region 
where the results are nearly independent of background 
assumed. We performed a number of tests to address 
the effects of the XMM-Newton background uncertain- 
ties on the deprojected temperature profiles. We did 
this by: 1) varying the blank-sky background normal- 
izations by ±10% (±20%) for the MOS (PN); 2) using 
the blank-sky backgrounds without renormalization; 3) 
using the MOSl±MOS2 spectra only; 4) using the PN 
spectra only; 5) adding an extra 0.2 keV MEKAL compo- 
nent (with normalization free and allowed to be negative 
in each annulus) as a model for any variation in the soft 
Galactic component. The test results are summarized in 
Figure [TU1 Our adopted deprojected temperature profile 
in this paper is shown in black diamonds. 

Within r < 180", all the deprojected temperatures 
are in very good agreement. All these tests show a sig- 



nificant jump in temperature near 240" and a flatten- 
ing of temperature beyond this jump. As expected, in- 
creasing the background normalization (blue diamonds) 
gave lower deprojected temperatures for the outer re- 
gions (r > 180"), while decreasing the background nor- 
malization (green diamonds) gave higher temperature. 
Using the blank-sky backgrounds without renormaliza- 
tion (red crosses) also gave higher temperatures (green 
diamonds). The variation in the temperatures using dif- 
ferent normalizations was at most ~ 1 keV in the outer 
region. 

The deprojected temperature profiles obtained from 
the MOSl±MOS2±PN spectra (black diamonds), the 
MOSl±MOS2 spectra (magenta triangles), and the PN 
spectra alone (cyan triangles) all agree reasonably well. 
Thus, the PN data appear to give fairly reliable spec- 
tral information, even though the uncertainty in the PN 
background normalization may be larger. 

Finally, varying the amount of soft Galactic back- 
ground (orange crosses) gave deprojected temperatures 
in good agreement with those using the renormalized 
blank sky backgrounds (black diamonds) except for the 
outermost annulus. The large deviation of the outermost 
data point may indicate that there is less of soft Galac- 
tic emission in the direction of Abell 2626 than in the 
blank-sky fields because the fitted normalization of the 
extra soft Galactic background component in the model 
is negative. The soft Galactic component should be uni- 
form on the scale of a cluster. However, a reduction in 
the soft Galactic background similar to that suggested by 
the last annulus was found to be completely inconsistent 
with the fits in the inner annuli. Hence, we concluded 
that the reduced soft component in the outermost an- 
nulus was not due to a reduction in the soft Galactic 
background towards the cluster. The reduced soft com- 
ponent in the outermost annulus may be the result of 
additional thermal structure in this region. 

To summarize, all our results are independent of back- 
ground used within ~ 180". For the outer regions, the 
temperature profile is only slightly affected by the choice 
of the background, but the differences are mostly well 
within the error bars. All the features described in the 
paper related to temperature measurements (e.g., the 
temperature jump, pressure, entropy profiles, etc.) are 
not affected. 

4.2. Central X-ray Spectrum and Cooling Rate 

We determined the radiative cooling rate by fitting 
the X-ray spectrum of the central regions of the clus- 
ter from both XMM-Newton and Chandra data. We 
extracted the spectrum from a circular region with a 
radius of 65", which is roughly the cooling radius de- 
termined in § 14.41 below. In XSPEC, we used the 
WABS*(MKCFLOW+MEKAL) model, which combines 
a cooling flow component (MKCFLOW) with an isother- 
mal model (MEKAL) to account for the gas in the outer 
regions. The absorption was fixed at the Galactic col- 
umn of 4.2 x 10 20 cm -2 as before. The upper tempera- 
ture (fcThigh) and abundance of the MKCFLOW compo- 
nent was set to be the same as the MEKAL component, 
which is expected if the gas cooled from the ambient 
ICM. Initially, the lower temperature of the MKCFLOW 
component is set to be the lower limit for the model 
(fcTi ow = 0.08 keV), as expected for the classical cooling 
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flow model in which the gas cools radiatively to very low 
temperatures. With this assumption, the best-fit cooling 
mass deposition rate is 2 ± 1 (3 ± 3) Mq yr _1 for XMM- 
Newton (Chandra) data in the 0.6-10.0 (0.5-7.0) keV 
energy band with a x 2 /dof = 1.18 (0.90). The fitted 
upper temperature and abundance were fc Thigh = 2.8 ± 
0.05 (3.3±0.2) keV and 0.52±0.03 (0.9±0.1) solar, respec- 
tively. Allowing kT\ ow to vary increased the mass deposi- 
tion rate to be 74±6 M Q yr" 1 for the XMM-Newton data, 
with fcT low = 1.5 ±0.1 keV, fcT high = 4.2^ keV, and an 

abundance of 0.45±q °2 solar, with a x 2 /dof = 1.13. The 
Chandra data gave very poor constraints on the mass de- 
position rate or lower temperature with a \ 2 /dof — 0.87 
if kT\ ow was allowed to vary. From the f-test, the proba- 
bility that the fit to the model with fcTi ow fixed is better 
than the fit with freeing fcTi ow is only 2 x 10 _11 (9 x 10~ 3 ) 
for the XMM-Newton (Chandra) data. For comparison, 
the mass deposition rate determined from Einstein data 
by surface brightn ess deprojection was ~ 53 M Q yr 



([White et al.l 1997|). As has been found in many other 



clusters (e.g. JPeterson et al"1l2003[ ). the spectral data for 
Abell 2626 can be fit with significant amounts of gas 
cooling, but only by a factor of ~3 in temperature. 

4.3. Density, Pressure and Entropy Profiles 

Assuming the cluster is spherically symmetric, which 
appears to be roughly true for Abell 2626, the radial 
surface brightness profile (Figure [6]) can be deprojected 
to give the X-ray emissivity e(r) and electron density 
n e (r) as a function of radius (Appendix |Aj . The elec- 
tron density profiles determined from XMM-Newton and 
Chandra data are shown in the upper panel of Figure [TT1 
Except for more fluctuations in the Chandra data, the 
two results for radii greater than ~ 7" are in good agree- 
ment. To avoid artifacts due to the assumption of zero 
emissivity outside the boundary during the deprojection 
(jTaki zawa eTaD[2003), we have removed the 2 outermost 
points in each of the profiles plotted in Figure [TT1 In the 
Chandra data, we also note that the two innermost data 
points are higher than the XMM-Newton results by a 
factor of ~ 1.5 — 2. Such a deviation in the innermost 
data points can also be seen when comparing the surface 
brightness profile of Chandra data to XMM-Newton data 
(not shown in this paper). The surface brightness profile, 
as well as other profiles, is centered at the location of the 
southwest nucleus of the cD galaxy. The central devia- 
tion between Chandra and XMM-Newton might be due to 
the X-ray contribution of a central AGN (§ 14. 7[) and the 
differences in the PSFs of XMM-Newton and Chandra. 
The smaller PSF of Chandra allows the central excess to 
be resolved more readily. We did not remove the central 
source during the analysis, since it is probably extended 
in the 0.3-10 keV band (see § 14. 7} . Also, including or 
excluding this region does not change any other conclu- 
sions in the paper. The central excess extends to 4" — 7" 
in radius, which is larger than the PSF of Chandra. This 
probably indicates that the excess is not due to an AGN 
alone. The northeast nucleus of the cD galaxy is located 
at about 4" away from the central AGN (Figure[3]). From 
the residual map (Figure [8]), we clearly see that there is 
an extended excess around the northeast nucleus. There 
is also excess emission northwest of the AGN, which can 
also be seen in the soft band raw image of Figure [3] 



The XMM-Newton data show that the electron density 
decreases from ~ 0.02 cm~ 3 at ~ 7" to ~ 0.0004 cm~ 3 at 
~ 400". The density profile is generally smooth, with a 
small jump around 60" seen in the Chandra data. How- 
ever, the jump is not visible with the XMM-Newton data. 
This radius corresponds to the location where the loga- 
rithmic slope of the surface brightness profile changes, or 
equivalently, to where the surface brightness goes from 
being dominated by the smaller to the larger of the two 
individual (3 models. We also noticed that the radio mini- 
halo is roughly confined within ~ 60". Interestingly, this 
radius is also roughly coincident with the cooling radius 
(§ 14.41 below). Within the cooling radius and the region 
of the radio mini-halo, radiative cooling and/or the ef- 
fects of the relativistic particles and magnetic fields in 
the radio plasma may affect the density profile of the 
X-ray gas. 

The gas density and temperatures profiles can be com- 
bined to give the gas pressure P = nkT, where n is the 
total number density in the gas (both ions and electrons). 
We assume a fully ionized plasma with half cosmic abun- 
dances of heavy elements which gives n ks 1.92 n e . We 

also derive the gas entropy parameter S = kT Jr?J^ 
(|Ponman et al.ll2003D . Since the spatial resolution of our 
density profile is much higher than the temperature pro- 
file, interpolation is used to determine the temperatures 
on the finer density grid by fitting a high order polyno- 
mial (lower panel of Figure [9]). We assigned tempera- 
ture errors of 5 (20)% at radii of smaller (larger) than 
250". Using a constant temperature model gives simi- 
lar results, except for more negative (unphysical) values 
for the total mass profile (§ 14.51 below). Thanks to the 
nearly constant temperature profile of Abell 2626, we be- 
lieve that the interpolation of the temperature profile is 
reasonably accurate. 

The pressure profile (middle panel of Figure [TT|) looks 
similar to the electron density profile, since the temper- 
ature is roughly constant. The feature at 60" becomes 
more obvious, as the slope of temperature is negative 
in this region. At around 240", there is a slope change 
in the pressure profile. This is associated to the tem- 
perature jump at that radius, which may be a merger 
signature. 

The entropy profile is shown in the lower panel of Fig- 
ure 111! There is a small bump in the entropy profile 
around 240" which is associated with the temperature 
jump and the change in slope of the pressure profile in the 
same region. The entropy may be boosted by a merger 
in that region. Excluding this entropy bump region, the 
logarithmic slope of the entropy profile in the outer re- 
gions (~ 100"-400") is 0.58±0.06, which is significantly 
small er than the universal scaling relation slope of ~ 1.1 
(e.g., iPonman et al. 1 120031 ). Note that there is a large 
scatter i n the entropy profiles of individual clusters (e.g., 
Figure l. lPonman et al.ll2003l ). and hence the smaller log- 
arithmic slope might not be unusual. There is a differ- 
ence in the central entropy measured by XMM-Newton 
and Chandra. This central entropy deviation comes from 
the difference in central gas density measurements, which 
is probably due to the combination of a central excess in 
X-ray emission and differences in the PSFs of Chandra 
and XMM-Newton. However, the spatial resolutions of 
the temperature measurements are not sufficient to de- 



X-ray Observations of Abell 2626 



9 



termine if the entropy differences are due to the tem- 
perature in the hot gas, or a central X-ray source. In 
the following, we chose to present results of the XMM- 
Newton data with which the spatial resolution in density 
profile measurement is poorer. The XMM-Newton data 
show that the entropy decreases from the outer value 
of ~ 360 keV cm 2 , eventually flattening to a value of 
~ 30 keV cm 2 . If there is no feedback mechanism, a pure 
radiative cooling model will result in a power law entropy 
profile, with central entro py lower than ~ 10 ke V cm 2 
within 10 kpc (Figure 1 of fVoit fc Donahue! l2005h . The 
central entropy pedestal value of ~ 30 keV cm 2 implies 
that there has to be heating in the central region, which 
can probably be explained by a central AGN. 

4.4. Cooling and Thermal Conduction Time Scales 

Three important processes, heating, cooling and heat 
conductio n, determine the thermal structure in the hot 
gas (e.g., ISarazin| [i986). While the heating process in 
the central region is probably due to the AGN, the other 
processes are related to the hot gas structure itself. Here 
we compare two important time scales for the hot gas 
- the cooling and thermal conduction time scales. For 
the temperature gradients seen in typical cooling flow 
clusters, conduction is likely to be an important process 
if it is not suppressed by the magnetic field. 

The cooling time was determined for each radius in 
the deprojected density profile using the values of the 
deprojected density and the deprojected temperature for 
both XMM-Newton and Chandra data. The rate of to- 
tal emission from the gas was determined using the same 
MEKAL model used to fit the spectra. The integrated, 
isobaric cooling time t coo i was determined. A brief out- 
line of the cooling time estimation procedure is given in 
Appendix [Bl 

We define the conduction time scale to be (|Sarazinl 
[19861) 

n e (T/\VT\n 

^cond — > \^) 

K 

where k is the ther mal conductivity for a hydrogen 
plasma (|Spitzerlll962f ). 

K = 4.6xl0 13 f^-Y /2 (—) 1 ergs s" 1 cm^ 1 K- 1 , 
V 10 8 KJ \ 40 J B 

(3) 

and A = 37.8 + ln[(T/10 8 K)(n e /10- 3 cm- 3 )- 1 / 2 ] is the 
Coulomb logarithm. 

Figure [12] shows the cooling and conduction time scales 
as a function of angular radius. The cooling time scale 
(plus symbols) decreases from ~ 50 Gyr at ~ 540", down 
to ~ 1 Gyr at < 15". Assuming the cluster has not 
changed much since z ~ 1, we define the cooling radius 
to be where t coo i is equal to the time since z = 1, which 
is 7.7 Gyr. The cooling radius, r coo i, is thus ~ 60" — 70" 
(~ 66 — 77 kpc). Interestingly, r coo i is of the same size as 
the mini-halo (see § 13. 4p . It is also located near the ra- 
dius where the logarithmic slope of the surface brightness 
profile changes, or equivalcntly, where the surface bright- 
ness profile goes from being dominated by the smaller to 
the larger of the two individual (3 models (§ 13 . 3[) . 

The conduction timescale depends on the temperature 
gradient, which is strongly affected by the uncertainties 



in the individual temperature points. The error estima- 
tion for the conduction timescale is complicated by the 
fact that the adjoining deprojected temperatures used to 
derive the conduction timescale are correlated. Hence, 
in calculating the uncertainties by error propagation, we 
included the full covariance matrix for the temperatures 
obtained from the deprojection using XSPEC. In the up- 
per panel of Figure [12l the crosses show the values of the 
conduction time derived using the individual tempera- 
ture values derived from annular XMM-Newton spectra 
(§ 14. 1[) . The conduction time derived from the Chandra 
data (lower panel of Figure !T2|) are in general shorter 
than those from the XMM-Newton data. This may be 
due to the much larger scattering of the Chandra depro- 
jected temperature profile. The error bars of the Chan- 
dra data are much larger than that of the XMM-Newton 
data. In order to reduce the effects of the uncertainties 
of the Chandra data, we also determine the conduction 
time scale in the inner region using the same interpo- 
lated temperature profile which was used to deproject 
the pressure/entropy profile (§ I4.3|) . The interpolated 
temperature profile was also used for the XMM-Newton 
data for comparison. These values are shown as dots 
in Figure [T^l the upper panel gives the values for the 
XMM-Newton data, and the lower panel is for the Chan- 
dra data. The Chandra conduction time scales are longer 
than the XMM-Newton values in the innermost regions 
due to the higher densities derived there from the Chan- 
dra data. 

Here we describe results obtained from XMM-Newton 
since the error bars are much smaller than those from 
Chandra. In the central regions of Abell 2626, the tem- 
perature gradient is very flat. The overall temperature 
variation is only about 0.4 keV within the cooling radius 
of ~65". Both the best-estimate and the lower bound on 
the conduction time are larger than the cooling time in 
the central regions of the cluster. The nearly isothermal 
temperature profile in Abell 2626 and resulting low con- 
duction imply that thermal conduction is not important 
compared to cooling anywhere the cooling rate is high. 
Hence, in Abell 2626 the thermal structure within the 
cooling radius of ~65" is probably determined by cool- 
ing and heating processes only. At a radius of ^240", the 
conduction time scale drops to ~ 0.3 Gyr only, which is 
significantly shorter than the look-back time to a red- 
shift z ~ 1. Even including the uncertainty in the con- 
duction timescale, the upper limit on the time scale is 
still quite short. This indicates that either the thermal 
structure (perhaps associated with a merger) happened 
within ~ 0.3 Gyr, or thermal conduction is suppressed 
at least by a factor of ^20. Except the feature at ~240", 
outside of ~65", we cannot determine whether thermal 
conduction is important or not due to the large uncer- 
tainties. 

4.5. Mass Profiles 

Figure Q2] shows the gas mass and total mass profiles 
of Abell 2626. The gas mass profile is given by 

M gas (<r)= f 47r(/) 2 d/p gas (r'), (4) 
Jo 

where p gas — [i e n e is the gas mass density, and ji e is the 
mass per electron determined from the gas abundances. 
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The total mass of the cluster can also be determined if 
we assume the hot gas is in hydrostatic equilibrium, and 
if the temperature of the ho t gas is know n. The total 
mass profile is given by (e.g.. ISarazinlll986T) : 

where [i is the mean mass per particle, and m p is the pro- 
ton mass. The temperature profile was interpolated from 
the XMM-Newton data (Figure [3 § |4~3|) . Since the mass 
within any given radius is determined by local gas prop- 
erties and their derivatives, the total mass can decrease 
or can even be negative. (This occurred at r = 57" for 
XMM-Newton and at a number of points for Chandra, 
and these points were omitted from the plot.) This can 
be due to measurement uncertainties in the data, which 
are compounded by the derivatives, or by the breakdown 
of simple hydrostatic equilibrium. The radio mini-halo 
is within ~ 60", and may be disturbing the hot gas or 
contributing additional pressure support in this region 
(§§ 13.1.21 & 13. 4[) . Thus, the mass profile at small radii 
may not be represented accurately by hydrostatic equi- 
librium. At radii of r > 60", the general trend is for 
M{< r) to incr ease with rad i us. 

Compared to IWhite et al.l ()1997l ) rescaled to our value 
of the Hubble constant, the gas masses we have ob- 
tained agree to within 2% and 7% at 321" and 489", 
respectively. The total mass agrees within 26% at 321". 
The difference in the total mass at this radius is large, 
but the values still agree within the uncertainties. The 
two values of the total mass do not agree at 489", 
probably due to the larger uncertainty in the temper- 
ature in the outermost annulus. The gas mass fraction 
(M gas /Mtot) increases from ~ 5% at 100" to ~ 10% at 
500". The glo bal trend is the same as is typical in re- 
laxed clusters (iDavid et alJl!995t lEttori fc Fabianlll999t 
lAllen et al.ll2002f ). 

4.6. Hardness Ratio Maps 

To understand the rough spectral properties of differ- 
ent regions, we have created hardness ratio maps from 
both XMM-Newton and Chandra data (Figures Q3] |T5l 
fc fT6]). For XMM-Newton data, the PN image was not 
included because of its different spectral response and 
because the chip gaps and bad columns produced cos- 
metic artifacts in the smoothed image. The hardness 
ratio maps were created by dividing the image in the 
hard band (2-10 keV) by the image in the soft band 
(0.3-2 keV). Each image was background-subtracted, 
exposure-corrected, and smoothed with the same scales 
as the 3cr adaptively smoothed image of the 0.3-10 keV 
band. We also tried smoothing all images with the same 
scales as the 5tr adaptively smoothed image of the 0.3- 
10 keV band, or with the same scales as the 3<r adaptively 
smoothed image of the hard band (2-10 keV), and almost 
all the features look the same as the hardness ratio maps 
given in this paper. While some features are not con- 
sistently seen in both the XMM-Newton and Chandra 
hardness ratio maps, indicating the possibility that they 
are statistic fluctuations in at least one of these maps, 
the common features which show up in both maps in 
the inner regions should be real (Figure [15]). The colors 
from black to blue to red to yellow to white represent the 



degree of hardness from soft to hard. Though the hard- 
ness ratio maps from XMM-Newton and Chandra cannot 
be directly compared due to different spectral responses, 
the relative hardness is still useful to identify interesting 
regions. 

4.6.1. Global Structure 

Figure [14] is the hardness ratio map of the whole 
Abell 2626 from XMM-Newton data, overlaid with green 
contours of excess emission from the XMM-Newton resid- 
ual map (left panel of Figure [7]). Part of the apparent 
inhomogeneity in hardness is just due to noise. Also, 
near the blue and dark boundary of the hardness ratio 
map, the features start to be dominated by the hard 
background. Although difficult to see due to the rather 
flat temperature profile, the overall picture is a soft core 
at the center. At larger radii, the hardness increases out- 
ward to a radius indicated around the central cyan (light 
blue) circle, and then the hardness decreases outward. 
The central cyan circle is at a radius of 240" where the 
jump in temperature was identified (§ I4.1|) . From the 
hardness ratio map, it seems that the harder gas is con- 
fined within the cyan circle. 

The white circles in Figure [M] are the two extended X- 
ray emission regions associated with the merging subclus- 
ters (§ 13.1. 1[) . The discontinuities in the green contours 
near the white circles are actually artifacts of chip gaps. 
Both of the extended X-ray emission regions appear to 
be harder than their surroundings, and each of them 
contains a softer point-like source. The soft source in 
the northeast extended X-ray emission region is located 
at the peak of the extended emission, and is probably 
associated with the galaxy 2MASX J23365372+2113322 
(Table [T]), while the soft source in the southwest region 
is not associated with any identified galaxy. The two 
soft point-like sources in the two extended emission re- 
gions disappeared when the images were smoothed with 
the same scales as the 5c adaptively smoothed image of 
the 0.3-10 keV band, or with the same scales as the 3a 
adaptively smoothed image of the hard band (2-10 keV). 

4.6.2. Central Region 

The details of the central region of the hardness ratio 
maps are shown in Figure [15] Contours of excess emis- 
sion are shown from the Chandra residual map (right 
panel of Figured]). While the XMM-Newton and Chan- 
dra hardness ratio maps do not agree completely, prob- 
ably due in part to noise and in part to the larger PSF 
of XMM-Newton, we can still identify some common fea- 
tures. 

One interesting feature is that the SO galaxy IC 5337 
shows a softer tail in the western direction. On the Chan- 
dra hardness ratio map (right panel of Figure [T5"]) , there is 
a narrow tail of soft emission starting from the southern 
corner of the IC 5337 tail and extending to the SWW. We 
also see a softer region with a similar scale in the XMM- 
Newton hardness ratio map (left panel of Figure [T5]) . 

The northern excess region ~ 30" above the cD galaxy 
(contours) appears softer, with an extension towards the 
northeast direction which can also be seen in Figure UM 
If the region is near hydrostatic equilibrium and the ma- 
terial is supported by ideal gas pressure, we would expect 
the brighter (and hence denser) regions to have softer 
spectra. This is roughly the case for the northern excess 
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region. However, for the southern excess at a similar ra- 
dius where the tongue is located (S I3.1.2j) . the brighter ex- 
cess region does not correspond to a softer region, which 
suggests a breakdown of hydrostatic equilibrium. In par- 
ticular, the X-ray excess tongue in the south direction of 
the cD galaxy (compared to Figure |HJ appears to over- 
lap both hard and soft regions extending in the same 
direction, making it hard to understand its origin. 

Between the main cluster and the southwest subclus- 
ter, the large extended region of excess emission near the 
southwest corner appears to have a harder spectrum in 
the Chandra data (right panel of Figure [15]), while it ap- 
pears to be soft in the XMM-Newton data (left panel of 
Figure [T5"|) . From the hardness ratio map, it is not clear 
whether the extended excess is a shock heated region or 
a cold front. 

In FigureHni we can see clearly that there are two hard 
point sources located at the position of the cD galaxy 
IC 5338 and the SO galaxy IC5337 (arrows). The inter- 
pretation of two AGNs will be discussed in § 14.71 Radio 
contours (|Gitti et all |2004) are overlaid on the Chan- 
dra hardness ratio map to show in detail the compli- 
cated thermal and radio structures in the central region. 
The southern radio bar appears to be in a harder re- 
gion, while the northern one appears softer. If the radio 
bars are indeed thin tubes along the plane of the sky, 
the X-ray emission along the line-of-sights may mainly 
be contributed by the cluster emission beyond the cen- 
tral region. In this case, the projected image or hardness 
ratio would not necessarily reflect the X-ray emissivitics 
or hardnesses of the radio bars, and hence the lack of 
symmetry or a clear correspondence with features in the 
X-ray images or hardness maps might be understandable. 

4.7. Possible X-ray AGNs in IC 5338 and IC 5337 

We searched for evidence of a central point source asso- 
ciated with the cD galaxy IC 5338, which may indicate 
the existence of an X-ray AGN. When wavedetect was 
run on the 0.3-10 keV and 0.3-2 keV band Chandra im- 
ages (§ I3.2p . each detected one central source located at 
the southwest and northeast nucleus of the cD galaxy, 
respectively. Both were possibly extended. We also ran 
wavedetect on the hard band (2-10 keV) image, but it did 
not detect a central point source. However, when we in- 
spected the hard band image, we noticed a possible hard 
point source. This source can be seen in the Chandra 
hardness ratio map (Figure fT6|) . We extracted the HST 
image of the center of IC 5338 from the archive. The 300 
s exposure was taken on 1999 June 26 using the HST 
WFPC2 camera with the F555W filter. Figure [3] shows 
clearly that this cD galaxy has two distinct nuclei, with 
the northeast nucleus being brighter optically. The two 
optical nuclei are located at the peaks of the soft X-ray 
image, with the radio core source in IC 5338 centered 
on the hard X-ray, southwest nucleus. The hard X-ray 
point source is only separated from the peak of the radio 
core by less than 1". Thus, we believe that this source 
is the AGN, and that the point-like hard X-ray source 
arises from this AGN. 

To study the X-ray properties of the IC 5338 AGN, 
Chandra data were used to do photometry. Due to the 
limited number of photons, we cannot fit the spectrum 
of the AGN. Hence, we estimated the properties of the 
cD AGN by assuming that it has a power law spectrum 



with a photon index T. We extracted the total pho- 
ton counts in a circular region with a 3 pixel radius in 
both hard (2-10 keV) and soft (0.3-2 keV) bands. In 
order to subtract the cluster emission as well as true 
background, we took background from a 3-pixel-wide 
annulus just outside the central region. The program 
PIMMS 6 was used to determine the photon index and the 
unabsorbed flux which gave the correct hardness ratio 
and photon counts. We assumed that absorbing column 
was given by the Galactic value Nh = 4.2 x 10 20 cm~ 2 
(|Dickev fc Lockmanlll990i ). The best fit photon index is 
r = 2.2±° 7 4 . The unabsorbed flux in the 0.3-10.0 keV 
band is Fx = 1.6^0 4 x 10~ 14 ergs cm -2 s _1 , which gives 
an X-ray luminosity of Lx = 1-2^03 x 10 41 ergs s _1 in 
the same energy band. 

We also searched for evidence for a point source as- 
sociated with an AGN in the center of the X-ray bright 
SO galaxy IC 5337. In this case, wavedetect did detect a 
point source in the Chandra hard band (2-10 keV) image 
at the center of IC 5337. This point-like, hard source is 
clearly seen in the hardness ratio map (Figure \W\i . The 
SO galaxy is also a radio source, and the X-ray point 
source is apparently associated with this central AGN. 
Again, the point source was too faint for spectral fitting. 
Using the same method described above, we determined 
the best-fit spectral index and flux assuming the Galactic 
hydrogen column density of 4.2 x 10 20 cm -2 . The best fit 
photon index is T — 0.2+Q g, the unabsorbed flux in the 

0.3-10.0 keV band is F x = 3.2l^xl0~ 14 ergs cm" 2 s~\ 
and the X-ray luminosity is Lx — 2.4l}'g x 10 41 ergs s -1 . 
This spectrum index seems to be low (hard) for a typ- 
ical AGN. One possibility is that there is excess in- 
ternal absorption associated with the SO galaxy AGN. 
To assess this, we fixed the photon index to T = 1.5. 
In this case, the best-fit hydrogen column density is 
N H = l-Stol x 1022 cm" 2 , the unabsorbed flux is 
F x = (3.0 ± 1.5) x 10~ 14 ergs cm~ 2 s -1 , and the X-ray 
luminosity is Lx — (2.3± 1.1) x 10 41 ergs s _1 . Somewhat 
counterintuitively, the extra absorption gives a lower un- 
absorbed flux because of the steeper spectral index as- 
sumed. 

5. KINEMATICS OF THE SO GALAXY IC 5337 

The bow-shock-like shape suggests that the SO galaxy 
is moving east directly towards the cluster center (right 
panel of Figurc[2]). The hardness ratio maps of the center 
of the cluster (Figures fT5lfc [T6|) show a region of soft (and 
presumably, cool) X-ray emission extending to the west 
of the SO galaxy. This geometry would be consistent with 
the SO galaxy falling into the cluster from the west; the 
cool tail of gas would be either interstellar or intragroup 
gas associated with this galaxy which is being stripped 
by ram pressure from the motion through the ICM of 
Abell 2626. The steepness of the intensity gradient in the 
cluster radio mini-halo on the western side compared to 
the eastern side (left panel of Figure |2J) may also indicate 
that it is compressed by the SO galaxy. 

The radio structure complicates the interpretation of 
the motion of the galaxy, since there is a component cen- 
tered on the SO galaxy, and a component located at a 
position of about 190° (measured from north to east) of 

6 http:/ /heasarc. gsfc.nasa.gov/Tools/w3pimms. html 
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the SO galaxy on the 1.5 GHz B-array map (Figures [2j 
fe[T6|). The 1.5 GHz C-array map shows three radio tails 
in the south, southwest and west directions. If the south 
component is a Narrow- Angle- Tail (NAT) behind the SO 
galaxy, then it is likely to be moving to the NNE (about 
10° measured from north to east) direction. On the other 
hand, the other two components (west and southwest 
tails) are more consistent with the bow-shocked structure 
in the X-ray image which indicates that the SO galaxy is 
moving to the east. 

The interpretation of the kinematics of IC 5337 is con- 
fused by an ambiguity in the literature concerning its 
radial velocity. There are two optical determinations 
of the radial velocity which give values which are con- 
sistent with the velocity of other galaxies in the cen- 
ter of Abell 2626, including the cD galaxy IC 5338. 
iHuchra et all (1 19991) gi ve a velocity of 16,485±37 km s -1 , 
and lFalco et alJ ([1999D find 16,562±60 km s _1 . These 
values are consistent wi th the 21 cm line determin ation 
of 16.485±40 k m s" 1 bvlGiovanelli fe Havnesl (fl993h . On 
the other hand. [Mohr et all (1996) l ist an optical veloc- 
ity of 18,903±39 km s _1 , and iKarachentsev fc Kopylovl 
(1981) list a velocity of 19,097 km s _1 with an un- 
certainty of the order of 100 km s _1 . The differ- 
ence between these velocities and the others is much 
greater than the quoted (or a ny likely) measure ment 
errors. The velocities given by iMohr et al.l (|1996l ) and 
IKarachentsev fc Kopylovl (|1981h would be inconsistent 
with IC 5337 being a part of the main central cluster 
Abell 2626. On the other hand, these higher velocities 
are consistent with IC 5337 being a member of the south- 
ern subcluster with a mean velocity of 19,164±138 km 
s _1 (|Mohr et al.l Il996h . Most of the other members of 
this subcluster lie considerably to the south. 

If IC 5337 is a member of the southern subcluster, its 
dynamics would probably follow somewhat the motion of 
the subcluster. We would expect IC 5337 would be in- 
falling from SSW to NNE (towards the position angle of 
30° direction measured from north to east). This would 
be inconsistent with our interpretation of the bow-shaped 
X-ray region to the east of IC 5337. Based on the X-ray 
image, we suggest that IC 5337 is falling directly into the 
main cluster from the west, and may not be associated 
with the southern subcluster. On the other hand, given 
that there are two velocity measurements co nsistent with 
the kinematics of the southwest subc luster ([Mohr et al.l 
11996. IKarachentsev fc Kopvlovi ri981). it seems unlikely 
that both are in error in the same way. One possible 
explanation is that there are two galaxies located along 
the line-of-sight toward IC 5337. Indeed, on the Chan- 
dra X-ray image, there are two point sources identified 
within the extended optical emission of IC 5337. One is 
located at the center of the SO galaxy with two west and 
southwest NATs, while the other is located at its south- 
ern edge with a south NAT. The southern X-ray source 
was removed during the analysis related to the SO galaxy 
in this paper. We have accessed the public available IR 7 , 
optical 8 and UV 9 surveys to search for a point source at 
the location of the southern X-ray point source. While 

7 2MASS: |http://irsa .ipac.caltech.edu/Missions / 2mass.html 

8 The STScI Digitize d Sky Survey: 
http: / /stdatu. stsci.edu/cgi-bin/dss_form 

a GALEX: http://www.galex.caltech.edu/ 



there is no obvious emission in the IR H and K band 
images, there seems to be extended emission associated 
with the southern X-ray source in the IR J band, optical 
band and the UV band images. However, one cannot 
easily determine if this emission is a separate galaxy or 
a feature in the disk of the SO galaxy. 

6. COOLING FLOW INTERACTION WITH THE CENTRAL 
RADIO SOURCE? 

The most unusual features of the radio image of 
Abell 2626 are the two elongated radio bars to the north 
and south of the radio core (Figures & [16]) . Their 
symmetric positions suggests they are radio lobes, but 
their elongated shapes are unusual. By comparison to 
radio lobes associated with other cooling core dominant 
radio galaxies, one might expect these to be regions of re- 
duced X- ray emission surrounded by bright rims ( "radio 
bubbles" . [Fabian et alj|200d : iBlanton et al.ll200lD . In ra- 
dio bubbles, the radio plasma has apparently displaced 
the X-ray gas. One would expect the level of reduc- 
tion in the X-ray surface brightness to depend on the 
extent of the bubbles along the line-of-sight. If they 
are highly oblate flattened regions with a large extent 
along the line-of-sight, the reduction in the X-ray sur- 
face brightness might be quite significant. If they are 
highly prolate tubes, then the reduction would be small. 
In fact, there is no obvious correlation between the two 
radio bars and the Chandra X-ray image, residual map 
or hardness ratio map in Figures & [TBI respectively. 
This may indicate that the radio bars are thin tubes, with 
a small extent along the line-of-sight, or that the radio 
plasma is mixe d with the X-ra y gas. In the numerical 
simulations by iRobinson et al.l (|2004T ). if the magnetic 
field of the bubble is weak, the bubble is very unstable 
and can be destroyed easily, leaving no obvious density 
contrast. The ratio of gas to magnetic pre ssure in the 
center of Abell 2626 is of the order of 100 (|Rizza et al.l 
12000. iGitti et all 120041 ), which is the sa me order as the 
weak magnetic field in the simulation by Robinso n et al.l 
(2004). However, it should be noted that the magnetic 
field in the radio bars could be much higher. One con- 
cern with this model is that it seems unlikely that such 
mixing would preserve the narrow radio bars which are 
observed. 

One possible interpretation of the arc-like X-ray fea- 
ture near the center of Abell 2626 (upper right panel of 
Figure [3]) is that the central cD galaxy is moving to the 
west relative to the local intracluster gas. This could ex- 
plain the bending of the X-ray emission near the center 
and the inner radio jets, as well as the elongation of the 
radio bars. However, in this case one might expect that 
the western edges of the two radio bars would indicate 
the positions where the radio jets have stopped. These 
positions are not mirror symmetric about the position of 
the radio core, as might be expected for jets perpendic- 
ular to the motion of the nucleus of the cD galaxy. Even 
if the jets are not perpendicular to the motion of the 
nucleus, one of the western edges of the two radio bars 
should have been behind the nucleus of the cD galaxy in 
this interpretation. 

There are two X-ray excess "tongues" from the cen- 
tral cD galaxy to the two radio bars, with the southern 
"tongue" appears to be stronger (upper right panel of 
Figure [3l & Figure [8]). A similar tongue has been seen in 



X-ray Observations of Abell 2626 



13 



Abell 133, where it was suggested that it might arise from 
a cluster merger, Kelvin-Helmholtz instabilities around 
the core, the buoyant uplift fro m a radio bubble, or a 
cooling wake (|Fuiita et al.l [2004") . Although subclusters 
may be merging with the outer regions of Abell 2626, 
there is no clear evidence for merger distortion in the 
central region. Even though the two cD nuclei may indi- 
cate a previous merger, the two nuclei are well inside the 
cD galaxy, which means that any associated merger may 
have happened a long time ago. Also, a merger might 
be expected to produce a tongue and radio lobe on only 
one side of the nucleus, as in Abell 133. For Abell 2626, 
the mean line-of-sight velocity of the galaxies is 16,533 
km s -1 , while that of t he central cD galaxy is 16,562 km 
s _1 (|Mohr et al.l 1199(1 ) . so the radial velocity of the cD 
galaxy relative to the cluster is low. Moreover, there is 
no obvious cold or shock front, which means that any 
motion in the plane of the sky is also slow. Therefore, 
there is no obvious velocity for the cD galaxy to induce 
Kelvin-Helmholtz instabilities. If the tongue is caused by 
a cooling wake, there would be only one. But we see two 
of them, though the northern one is much weaker than 
the southern one. 

Eliminating the explanations by a cluster merger, 
Kelvin-Helmholtz instabilities and the cooling wake, the 
tongue seems to be best explained by the buoyant up- 
lift from a radio bubble. In fact, a recent simulation 
indicates that when a bubble rises in the intraclustcr 
medium, large rolls extend out and converge under the 
bubble. This results in an upwelling and compression of 
mate rial under the wake of the bubble (Robinso n et al.l 
2004). The morphology of the southern tongue shown 
in Figure [H looks ver y similar to the simulated result 
(jRobinson et al.ll2004h if we assume that the radio bar is 
the head of the plume. However, if the southern "tongue" 
is caused by uplifted cool gas, why is the northern X- 
ray tongue associated with the northern radio bar much 
weaker, given the two radio bars are similar in appear- 
ance? While there is no radio flux measurement to the ra- 
dio bars in the literature, the stronger of the radio emis- 
sion from the southern bar can be noted on t he VLA 330 
MHz B +DnC-array radio map of Figure 4 in lGitti et al.l 
(|2004h . This may explain the stronger uplift of the ther- 
mal plasma in the south. 

Jet precession might also provide a natural explanation 
of the structure of the central region. The southwest 
cD nucleus is an AGN (see § I4.7j) . and may be eject- 
ing two jets towards the north and south. We assume 
the two jets are precessing about an axis which is nearly 
perpendicular to the line-of-sight along the north-south 
direction, so that the jets sweep out two conic surfaces 
to the north and south. If the two jets are stopped at 
approximately equal radii from the AGN (at a "working 
surface" ) , and if radio emission is produced by particles 
accelerated when the jets are stopped, then the radio 
bars might be produced as the jets precessed. If the 
jets are narrow, then the radio emission at the work- 
ing surface may be narrow, explaining the shape of the 
radio bars and also the lack of anticorrelation between 
the two radio bars and the X-ray emission. The preces- 
sion of the jets might be due to the gravitational effects 
of the second cD nucleus to the north, and this nucleus 
might have disrupted or weakened the northern X-ray 
tongue. Alternatively, the dominance of the southern X- 



ray tongue may simply be due to the southern jet being 
more powerful. There is also evidence of jet p recession 
in other cooling core clusters (jGitti et al.ll2006h . Jet pre- 
cession has also been suggested to be a requirement for 
the AGN feedback mechanism to so lve the cooling flow 
problem (jVernaleo fc R eynolds 200(|). 

Abell 2626 is one of the clusters associated with a pos- 
sible radio mini-halo, an extended diffuse radio source 
in the cooling co re region of the intracluster medium 
(jGitti et al.ll2004h . In the cooling core region, the mag- 
netic field strength is expected to be high due to com- 
pression. Consequently the radiative cooling time of 
the relativistic electrons is expected to be as short as 
~ 10 7 — 10 8 yr, which is much shorter tha n the time 
neede d to travel across the extended region (jGitti et al.l 
l2004h . In order to explain the existence of the radio mini- 
halo, an in-situ re-acceleration mechanism is needed. 
The steep radio spectrum of the mini-halo (a ~ 2.4, 
iGitti et al.ll200l may be due to confinement by the X- 
ray gas and radiative losses by the relativistic electrons. 

While there is no very obvious correlation between the 
radio emission from the mini-halo and the X-ray emission 
from the X-ray images (Figure [2|), there is some evidence 
of interaction between the radio mini-halo and the sur- 
rounding intracluster medium in the residual maps (Fig- 
ure [8]). In general, the diamond-shaped outline of the 
extended radio emission (the radio "mini-halo") corre- 
sponds to a region of reduced X-ray emission, and there 
are excess X-ray regions NNE and SSW of the edge of 
the mini-halo where the radio surface brightness is drop- 
ping rapidly. The mini-halo is roughly confined within 
60". The distance from the center to the NE or SW cor- 
ners of the diamond is about 50"; while to the NW or 
SE corners is about 70 ". As noted prev iously (§ 13. 4| . 
in the interpretation of Gitti et al. (2004) the two radio 
bars are distinct from the diffuse radio emission, and in 
fact after their subtraction the morphology of the mini- 
halo become roughly circular (see Figure 6 of lGitti et al.1 
2004), though we can still see some excess radio emis- 
sion in the NW and SE corners. Interestingly, there is 
also a small density jump in the X-ray gas around 60" 
(§ 14.31 above). The size of the mini-halo is also coinci- 
dent with the cooling radius (§ l4.4l above). This supports 
the re-acceleratiq n model for the origin of the mini-halo 
(jGitti et al.ll2004D . though it could just be a coincidence 
since the definition of the cooling radius is somewhat ar- 
bitrary. 

An excess in the residual map corresponds to a region 
where the gas density is higher than in the symmetric 
model (Figure [8]). The shape of the mini- halo, if it is re- 
ally elongated in shape (after proper subtraction of the 
radio bars), may be due to variations in the X-ray gas 
density. The radio source may be more strongly confined 
by the gas in the NE and SW directions. The points on 
the diamond-shaped radio mini-halo to the NW and SE 
may indicate regions of lower gas density where the ra- 
dio plasma has been able to expand more. On the west- 
ern side of the radio mini-halo, the intensity gradient is 
steeper than that of the eastern side. This may be due 
to gas compression by the SO galaxy IC 5337 falling to- 
wards the cD galaxy IC 5338. If this is the case, then 
IC 5337 should be very close to the central cooling core, 
rather than being part of the southwest subcluster. This 
reinforces both our dynamical picture for the SO galaxy 
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and the interpretation of the mini-halo geometry by con- 
finement. The X-ray deficient regions in the NW and 
SE directions can also be seen on XMM-Newton residual 
map on a larger scale (dark regions of the left panel of 
Figure [7|) . The excess of intracluster gas density to the 
NE and SW might be due to mergers along the axis of 
the local large scale structure, as discussed in § 13.1.11 or 
it may just represent the elongation of the equilibrium 
distribution of the gas in this direction. The lower den- 
sity regions in the NE and SW directions can also be 
regions inflated by previous radio bubbles. Although the 
symmetric radio bars may be due to the jet precession, 
it is also possible that the denser gas to the NE and SW 
might explain the shape of the two radio bars; the radio 
bars might have been compressed radially by a collision 
with this gas, and become extended to the NW and SE 
by the smaller ICM density gradients in those directions. 

7. CONCLUSIONS 

We have identified two extended X-ray emitting re- 
gions about 7' northeast and southwest from the center 
of Abell 2626, with the latter one probably associated 
with a subcluster identified to be falling into the main 
cluster (|Mohr et all H99l iMohr fc Wegnerl 119971 ). We 
argue that the infalling velocity of the southwest sub- 
cluster is mainly along the line-of-sight. Abell 2626 is 
known to be associated with the P erseus-Pegasus su- 
percluster (iBatuski fc Burnslll985al lbl; IZucca et alj|1993t 
lEinasto et al.ll2001l ). and it is likely that the main cluster 
is accreting sub-clusters and groups from this large-scale 
structure (§§ 13.1.1113.41 fc l4.6.2p . From the hardness ratio 
map, both of the extended X-ray emitting regions appear 
to be harder than their surroundings, and each of them 
contains a softer point-like source. The northeast point- 
like source may be associated with the galaxy 2MASX 
J23365372+2113322. 

Based on the bow shape of the X-ray emission and the 
soft X-ray tail identified by the hardness ratio map, we 
argue that the SO galaxy IC 5337 is falling towards the 
center of the main cluster from the west. The steeper 
intensity gradient in the cluster radio mini-halo on the 
western side compared to the eastern side may also in- 
dicate that the radio plasma is compressed by the SO 
galaxy. This reinforces our interpretation that the SO 
galaxy should be close to the cooling core region instead 
of located inside the southwest subcluster. On the other 
hand, if the SO galaxy were associated with the subcluster 
in the southwest direction, this motion would be hard to 
understand. While there are two velocity measurements 
of IC 5337 consistent with the southwest subcluster and 
three measurements consistent with the main cluster, one 
possible explanation is that there may be two galaxies 
along the line-of-sight of IC 5337. 

We have identified hard X-ray point sources due to 
AGNs associated with the cD galaxy IC 5338 and the SO 
galaxy IC 5337, both of which are radio sources. The 
cD galaxy AGN is located on the southwest nucleus of 
the cD galaxy. The SO galaxy AGN shows evidence for 
excess internal absorption. 

We have used a double beta model to fit the radial sur- 
face brightness profile and generated residual images to 
understand the dynamics of the intracluster gas (§ I3.4|) . 
Most of the residual structures in the inner region occur 
on a scale similar to that of the radio sources in the cen- 



ter. One major residual in the inner cluster which is not 
on the scale of the radio emission is the extended excess 
to the southwest. It is located between the cD galaxies 
and the southwest subcluster, and may be the result of 
a merger. 

The gas temperature in Abell 2626 increases from a 
central value of ~ 2.5 keV at ~ 10 kpc to a maximum 
of ~ 3.5 keV at around 230 kpc, and then decreases to 
~ 1.5 keV at ~ 600 kpc. At a radius of ~ 260 kpc, 
there is a significant temperature jump seen in the fits 
to the XMM-Newton spectra. The jump may be asso- 
ciated with a previous or an ongoing merger. We used 
the X-ray spectra to determine the cooling rate within 
a projected radius of 72 kpc. The amount of gas which 
is cooling down to very low temperatures is less than a 
few Mq yr" 1 . A much higher mass deposition rate of 
~ 74 M Q yr" 1 is obtained from the XMM-Newton spec- 
tra if the hot gas only cools by about a factor of three. 

From the surface brightness and the temperature pro- 
files, we have derived the density, pressure, and entropy 
profiles of Abell 2626. At a radius of ~ 70 kpc, there 
is a slope discontinuity in the density and pressure pro- 
files. The outer entropy profile has a logarithmic slope of 
~ 0.58, which is lower than that of the universal scaling 
relation of ~ 1.1. A heating source is probably required 
to maintain the central entropy of ~ 30 keV cm 2 . At a 
radius of ~ 260 kpc, there is a small bump in the entropy 
profile which may be associated with a merger. 

Due to the nearly isothermal structure of Abell 2626 
within the cooling radius, conduction is not important 
compared to cooling in the inner region. The cooling ra- 
dius is ^70 kpc, where the logarithmic slope of the sur- 
face brightness profile changes. At a radius of ~ 260 kpc, 
the conduction time scale is significantly shorter than the 
look-back time to a redshift of z ~ 1, which indicates that 
either the thermal structure (perhaps associated with 
a merger) happened recently or thermal conduction is 
highly suppressed. We have also derived the total mass 
profile under the assumption of hydrostatic equilibrium, 
and compared the gas mass fraction with other relaxed 
clusters. The global trend is typical for relaxed clusters. 

We have investigated the cooling flow interaction with 
the central radio source. The two elongated radio bars to 
the north and south of the center of the cD galaxy and 
the lack of obvious correlation between the two radio 
bars and any structures in the Chandra X-ray image, 
residual image, or hardness ratio map may indicate that 
the radio bars are thin tubes parallel to the plane of the 
sky. Another possibility is that the radio plasma is mixed 
with the X-ray gas, rather than displacing it. 

The central radio structure may be due to jet preces- 
sion. The southwest cD nucleus is an AGN, and we sug- 
gest it has two precessing jets propagating towards the 
north and south, producing the two symmetric narrow 
radio bars. The precession might be caused by the grav- 
itational influence of the second cD nucleus to the north. 
The two jets may have uplifted cool gas from the cen- 
tral region, producing the X-ray excess "tongues" run- 
ning from the center of the cD galaxy to the radio bars. 
The northern tongue appears to be weaker, and it may 
be disrupted by the northern nucleus of the cD galaxy. 
Alternatively, the strength of the southern X-ray tongue 
may simply be due to the southern jet being stronger. Jet 
precession may also be found in other cooling core clus- 
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ters (jGitti et al.ll2006[ ). It may help to solve the AGN 
feedback problem found in high resolution simulations 
(|Vernaleo fc Reynolds! [2006) . Alternatively, the distor- 
tions of the radio source might be due to ICM motion, 
particularly rotation, near the center of the cluster. 

The radio mini-halo has a diamond shape, elongated 
in the northwest and southeast direction. The diamond- 
shaped radio mini-halo seems to be compressed from 
the northeast and southwest direction, which is proba- 
bly consistent with the merging scenario for Abell 2626. 
With this compression geometry, the radio mini-halo 
might be leaking out in the northwest and southeast di- 
rection (Figured]). The western side of the mini- halo may 
also be compressed by the infalling SO galaxy IC 5337. 
Interestingly, the size of the mini-halo is roughly the 
same as the cooling radius, and is located at where 
the logarithmic slope of the surface brightness profile 
changes. The agreement between the size of the ra- 
dio mini-halo and the cooling radius is consistent with 
the re-acceleration model for the origin of the mini-halo 
(jGitti et al.ll2004D . although it might just be a coinci- 
dence since the definition of the cooling radius is some- 
what arbitrary. Independent of the origin of the radio 



mini-halo, the change in the X-ray surface brightness 
slope at the outer edge of the radio mini-halo may in- 
dicate that the radio plasma is affecting the thermal 
plasma. 
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APPENDIX 
DENSITY DEPROJECTION 

We basically followed iKriss et al.l (|1983f ) to deproject the gas distributions from surface brightness profile. The 
procedure is briefly outlined here. 

Assuming spherical symmetric gas distribution, the observed X-ray surface brightness (in ergs cm~ 2 sec -1 Hz -1 
sr _1 ) at an angle 9 (corresponding to a projected physical radius b) from the center of the cluster is: 

1 f°° e v ,{r)dr 2 



I° hB (9) 



4tt(1 + z) 3 J b 2 {r 2 -6 2 ) 1 /2 



(Al) 



where e„ is the X-ray emissivity of the gas integrated over all angles, v is the observed frequency, v' = (1 + z)v is the 
emitted frequency, z is the redshift, and r is the physical distance from the cluster center. 

We divided the image into n annuli bounded by the angular radii 9q < 9\ < ■ ■ ■ < 9i < 9i + \ <•••<#„ and centered 
on the peak of the cluster emission. In this paper, we take 9q = 0. Noting that b = 9dA, where dA is the angular 
diameter distance, the flux from the annulus from — > 9 i+ i is given by: 



F 



obs / 



J i+l, 



I° hs {9) 2tt 9 d9 



d A J '•> 



I° hs (b) db l 



(A2) 



Putting equation (jAljl into (IA2j) . discrctizing the integrals into summations, and assuming a constant emissivity 
(e v 'j) within each radius (bj — > foj+i), the observed flux can be expressed as: 



F 



obs / 



E ^M + i-%? ,2 -^-%?' 2 



j=i+i 



-(^ +1 -^i) 3/2 + (^-^i) 3/2 ]} ; 



or in matrix form: 



F° bs (9 l 



(A3) 



(A4) 



where Aij contains the geometrical factors in equation (|A3|) . The matrix is triangular, and equation (|A4|) can 
easily be inverted to give the emissivities at each radius. Once the emissivity is known, it can then be converted to 
the electron density, n e , by the relation: 

t v = n 2 A„ , (A5) 

where A„ is the emissivity function which depends on the temperature and the abundances. Note that A^ can be 
derived from the models used to fit the X-ray spectra. 
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COOLING TIME ESTIMATION 

The cooling time of a region can be estimated from the normalization of the MEKAL model fitted by XSPEC, 
together with the fitted temperature, the bolometric luminosity estimated by XSPEC, and the deprojected density. 
The normalization of the MEKAL model in XSPEC is given by: 

1CT 14 f 

K= H(iT? j " gMycm " ' (B1) 

where d,A and z are the angular diameter distance and the redshift to the source, respectively, n e and n# are the 
electron and hydrogen densities, respectively and V is the volume of the region. It is easy to show that the enthalpy 
H of the region is: 

H=-PV 
2 

2.5 x 2.33 x 4tt4(1 + z) 2 W u K , _ 
w ^ '- fcTergs. (B2) 

n e 

The instantaneous isobaric cooling time is then given by: 

Col = 7 — > ( B3 ) 

-^bol 

where Lboi is the bolometric luminosity of the region which is given by XSPEC. The bolometric luminosity can be 
written as £t>oi = MT)n1V, where the emissivity function A(T) depends on the abundances and is taken from the 
same XSPEC model which best fits the spectrum. Then, the integrated cooling time is 

A(T) [ T T'dT _ t* cool 



icooi-Coi T 2 j Q A(T') ~ 2° ' (B4) 
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TABLE 1 

XMM-Newton X-ray Point Sources with NED Identifications 



X-ray RA 


X-ray Dec 


NED object 


Type" 


NED RA 


NED Dec 


Redshift 


Separation 1 ' 


23:36:04.8 


+21 


:06:11.2 


2MASX J23360455+2106127 


G/IrS 


23:36:04.6 


+21 


:06:13 


0.066130 


3" 


23:36:09.6 


+21 


:20:25.1 


MG3 J233615+2120 


RadioS 


23:36:09.7 


+21 


:20:25 




1" 


23:36:15.8 


+21 


:06:11.5 


Abell 2626:[SPS89] 08 


VisS/RadioS 


23:36:16.2 


+21 


:06:21 




11" 


23:36:24.5 


+21 


:09:03.6 


IC 5337 


G/IrS/RadioS 


23:36:25.1 


+21 


:09:02 


0.054988 


9" 


23:36:30.5 


+21 


:08:47.0 


IC 5338 


G/IrS/RadioS 


23:36:30.6 


+21 


:08:50 


0.054900 


3" 


23:36:53.8 


+21 


:13:32.9 


2MASX J23365372+2113322 


G/IrS 


23:36:53.7 


+21 


:13:32 


0.056100 


2" 


23:36:54.2 


+21 


:15:31.7 


2MASX J23365415+2115302 


G/IrS 


23:36:54.2 


+21 


:15:29 


0.054200 


3" 


23:36:57.1 


+21 


:14:05.6 


2MASX J23365722+2114032 


G/IrS 


23:36:57.2 


+21 


;14:03 


0.038040 


3" 


23:37:32.4 


+21 


:09:07.6 


2MASX J23373222+2109020 


G/IrS 


23:37:32.2 


+21 


:09:02 


0.184400 


6" 



a The source type as classified by NED: G— Galaxies; IrS— Infrared Source; RadioS— Radio Source; VisS— Visual Sources. b Separation 
between the X-ray point source and source identified by NED. 



TABLE 2 
Chandra X-ray Point Sources 



RA Dec 



23 


:36 


:15 


.82 


+21: 


:08: 


:55 


.4 


23 


:36 


:21 


.31 


+21: 


13 


:29 


,7 


23 


:36 


:21 


.41 


+21: 


:11: 


35 


.4 


23 


:36 


:24 


.61 


+21: 


:08: 


.47. 


.7 


23 


:36 


:24 


.81 


+21: 


:10 


:16. 


,7 


23 


:36 


:25 


.04 


+21: 


:09 


03 


.0 


23 


:36 


:25 


.60 


+21: 


:13: 


:07 


,7 


23 


:36 


:30 


.52 


+21: 


:08 


:47. 





23 


:36 


:33 


.30 


+21: 


:08: 


:36. 


,7 


23 


:36 


:36 


.98 


+21: 


:10: 


:51. 


.8 


23 


:36 


:37 


.46 


+21: 


:06: 


:21 


.0 


23 


:36 


:38 


.43 


+21: 


:08 


:26. 


1 


23 


:36 


:15 


.60 


+21: 


:07: 


:49 


3 


23 


:36 


:46 


.84 


+21: 


:08: 


:07. 


.8 
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Fig. 1. — Background-subtracted, exposure-corrected, adaptively smoothed mosaic of the XMM-Newton EPIC MOS1 and MOS2 images 
in the 0.3-10 keV band. The image was smoothed to a signal-to-noise ratio of 3. The units for the color scale are counts per binned pixel, 
where each binned pixel has a size of 4'.'1 X 4'.'1. The arrows on the left and right indicate the cD galaxy IC 5338 and the SO galaxy IC 5337, 
respectively. Two outer extended emission regions are indicated with white circles. It should be noted that there is a chip gap at about 
5.'5 from the center which makes the northeast extended X-ray source appear to be more distinct than it actually is. 
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Fig. 2. — Left panel: Background subtracted, exposure-corrected, adaptively smoothed Chandra image of the central 3f5x2f4 region 
of the cluster in the 0.3—10 keV band. The image was smoothed to a signal-to-noise ratio of 3. The Chandra image shows the details of 
the inner region of the cooling core. The white solid lines are 1.5 GHz C-array radio contours showing the mi ni-halo. Th e green solid 
lines are 1.5 GHz B-array radio contours showing the radio bar structures. Both radio contours were taken from lGitti et alj 0001) . The 
arrows on the left and right indicate the cD galaxy IC 5338 and the SO galaxy IC 5337, respectively. Right panel: Background subtracted, 
exposure-corrected, adaptively smoothed Chandra image centered on the SO galaxy IC 5337. The image was smoothed to a signal-to-noise 
ratio of 4, with an intensity scale chosen to show the bow-shock-like shape. The southern source (red ellipse) was removed and replaced by 
the average surrounding X-ra y intensity before adaptively smoothing. The green contours are from the VLA 1.5 GHz B-array radio image 
taken from Gitti et al. (2004), with the lowest contour level being a signal-to-noise ratio of 3. 
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Fig. 3. — Images of the two nuclei of the cD galaxy IC 5338 in Abell 2626. All images show the same field of view of 30" X 30". Upper 
left panel: Optical image from HST archive. The image was taken with WFPC2 using the F555W filter. The two green circles, included in 
the other panels, are centered at the two optical nuclei observed by HST. Upper right panel: Background-subtracted, exposure-corrected, 
adaptively smoothed Chandra image in 0.3-10 keV band. The image was smoothed to a signal-to-noise ratio of 3. The color represents 
the X-ray intensity from high (white yellow) to low (dark blue). Lower left panel: Raw Chandra image in the soft (0.3-2 keV) band. 
Two intensity peaks can be identified. Lower right panel: Raw Chandra image in the hard (2-10 keV) band. Only the southwest nucleus 
corresponds to the peak in the hard band. 
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Fig. 4. — Optical image from the Digital Sky Survey showing the 9 XMM- Newton X-ray sources with identifications in NED (Table[T]l. 
The red circles are centered on the X-ray positions determined with XMM-Newton. There are green circles centered on the optical 
positions from NED, but they overlap so completely with the red circles that they are barely visible. NED objects are considered possible 
identifications if they are within 20" of an X-ray source. In fact all optical counterparts identified are within 11". 
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Fig. 5. — Optical image from the Digital Sky Survey showing the positions (red circles) of the 14 X-ray point sources detected with 
Chandra (Table 0. 
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Fig. 6. — Upper panel: Azimuthally averaged X-ray surface brightness profile from the XMM-Newton MOS1 and MOS2 observations 
of Abell 2626. The squares are the X-ray data, while the solid line is the fit to a double f3 model. The tiny error bars within the squares 
on the surface brightness are too small to be easily seen. Lower panel: The residuals to the best fit double beta model, in unit of <r. The 
horizontal error bars show the widths of the bins used to accumulate the counts. The vertical error bars are at the 1 a level. 
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Fig. 7. — Left panel: Residual map from the XMM-Newton MOS1 and MOS2 observations of Abell 2626, generated from a background- 
subtracted image. Yellow represents the largest excess, while black is the greatest deficit. Right panel: Residual map from the Chandra 
observation of Abell 2626, generated from a background-subtracted image. Yellow represents the largest excess, while black is the greatest 
deficit. There is an excess coincident with the center of the central cD galaxy, which is partly surrounded by a trough of reduced X-ray 
surface brightness at a radius of about 15". There are two localized regions of excess south and north of the central peak, which are 
indicated with green circles. There is an extended deficit about 50" west of the center, which is indicated with a white polygon. An 
extended region of excess emission is seen about 110" SW of the center, which is indicated with a green polygon. 
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Fig. 8. — Residual map from Chandra observation of the central 3.'5 X 2.'4 region of Abell 2626, generated from a background-subtracted 
image. The color scale from yellow to red to blue to black represents the residuals, running from excess to deficit. The white solid lines 
are 1.5 GHz C-array radio contours showing the mini-halo. The gre en solid lines are 1.5 GHz B-array radio contours showing the radio bar 
structures. Both radio contours were taken from lGitti eTaTI 112001 ). The arrows on the left and right indicate the cD galaxy IC 5338 and 
the SO galaxy IC 5337, respectively. 
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Fig. 9. — Upper Panel: Profile of the azimuthally averaged, projected temperature T of Abell 2626. The XMM-Newton (Chandra) 
data are shown as solid (open) circles. Middle Panel: Profile of the azimuthally averaged, deprojected temperature T of Abell 2626. The 
XMM-Newton (Chandra) data are shown as solid (open) circles. The black crosses are the projected temperatures from the XMM-Newton 
data in the upper panel for comparison (without error bars for clarity). In both the upper and middle panels, the error bars are at the 
90% confidence level. Lower Panel: Crosses are the temperatures calculated from interpolation of the original XMM-Newton deprojected 
te mperatur e p rofil e; this interpolated profile was used for the pressure/entropy/time scales deprojection and the mass profile analysis in 
§§ 14-31 14.41 & 14.51 The interpolated temperature profile for the Chandra data analysis is similar but not shown for clarity. The errors 
are assigned as 5 (20)% at radii of smaller (larger) than 250" to include possible systematic errors. The solid circles are the original 
XMM-Newton measured deprojected temperatures as shown in the middle panel, without error bars for clarity. 
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Fig. 10. — Profile of the azimuthally averaged, deprojected temperature T of Abell 2626 with different XMM-Newton backgrounds. The 
black diamonds give our adopted XMM-Newton temperature profile (the black solid circles in the second panel of Figure |9). The other 
temperature profiles are generated by: la) (blue diamonds) increasing the blank-sky background normalizations by 10% (20%) for MOS 
(PN); lb) (green diamonds) decreasing the blank-sky background normalizations by 10% (20%) for MOS (PN); 2) (red crosses) using the 
blank-sky backgrounds without renormalization (no error bars for clarity); 3) (magenta triangles) using the MOS1+MOS2 spectra only 
(data shifted to the right for clarity); 4) (cyan triangles) using the PN spectra only (data shifted to the left for clarity); 5) (orange crosses) 
adding an extra 0.2 keV MEKAL component to represent variations in the soft Galactic component (normalization free and allowed to be 
negative in each annulus, data shifted to the left for clarity). The error bars are at the 90% confidence level. 
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Fig. 11. — Upper Panel: Azimuthally averaged electron density profile of Abell 2626. The XMM-Newton (Chandra) data are shown 
as black solid (red dashed) symbols. For clarity, we have also multiplied the Chandra data and error bars by a factor of 2 (green solid 
symbols). The error bars are at the 1 a level. Middle Panel: Azimuthally averaged pressure profile, with the same notation as the upper 
panel. Lower Panel: Azimuthally averaged entropy profile, with the same notation as the upper panel. 
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Fig. 12. — Upper Panel: Cooling and thermal conduction time scales in Abell 2626 obtained from the XMM-Newton data. The cooling 
times are the plus symbols with horizontal bars showing the annular regions sampled. The conduction time scales derived directly from 
individual data points are the crosses with larger vertical error bars. The dots are conduction time scales calculated by interpolating the 
temperature profile of the cluster (see text in § 14.41 1 . The radial regions for the conduction time are not given for clarity. The error bars 
are at the 1 a level. Lower Panel: Same as the upper panel, but with the time scales calculated from the Chandra data. 
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Fig. 13. — Upper panel is the interior mass profile [M(< r)], of Abell 2626 obtained from XMM-Newton data, while the lower panel is 
from Chandra data. The total mass is shown with diamond symbols, while the gas mass points are circles. The error bars on the masses 
are calculated using propagation of errors at 1 a level. The error bars on the gas mass are plotted, but are smaller than the circle symbols. 




Fig. 14. — Hardness ratio (HR) map of Abell 2626 from the XMM-Newton observation. Background was subtracted from the images 
combined to calculate the hardness ratio. The color scale ranges from black (HR ~ 0) to blue (HR ~ 0.1) to red (HR ~ 0.3) to yellow 
(HR ~ 0.5) to white (iTi? ~ 2.5), showing increasingly hard emission. The green solid lines are contours from XMM-Newton residual map 
(left panel of Figure 01 showing positive residuals (excess emission). The white circles are the regions of extended X-ray emissions to the 
northeast and southwest. The cyan circle has a radius of 240" corresponding to the jump in the temperature profile. 
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Fig. 15. — Le ft p anel: Central portion (4.'6 X 4'.7) of the XMM-Newton hardness ratio map (Figure |14> . The color scale is similar 
to that in Figure IT4l but with the hardness ratio extremes being ~ 0.2 and 0.35 (dark blue to white yellow). The green solid lines are 
contours from the Chandra residual map (right panel of Figure [3 showing the excess emission. Right panel: Hardness ratio map from the 
Chandra data, for the same regions and with the same contours as the left panel. Background was subtracted from the images combined 
to calculate the hardness ratio. The color scale is similar to that in Figure [14] but with the hardness ratios in the range from ~ 0.1 to 
1.5 (dark blue to white yellow). The cD galaxy IC 5338 and the SO galaxy IC 5337 are located at (23:36:31, +21:08:47) and (23:36:25, 
+21:09:03), respectively. 
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Fig. 16. — Central region (3.'5 X 2f4) of the hardness ratio map from the Chandra observation. Background was sub tracted from the 
images combined to calculate the hardness ratio. The color scale is roughly the same as in the right panel of Figure 1151 The white solid 
lines are 1.5 GHz C-array radio contours showing the mini-halo. The gr een solid lines are 1.5 GHz B-array radio contours showing the radio 
bar structures. Both radio contours were taken from lGitti eTaTl 12004) . The arrows on the left and right indicate the cD galaxy IC 5338 
and the SO galaxy IC 5337, respectively. 



